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The aim of this chapter is to summarise the basic arguments and the intuition bolstering the 
RFOT picture for glasses, based on a finite dimensional extension of mean-field models with an 
exponentially large number of metastable states. We review the pros and cons that support or 
undermine the theory, and the directions, both theoretical and experimental, where progress is 
needed to ascertain the status of RFOT. We elaborate in particular on the notions of mosaic state 
and point-to-set correlations, and insist on the importance of fluctuations in finite dimensions, that 
significantly blur the expected cross-over between a Mode-Coupling like regime and the mosaic, 
activated regime. We discuss in detail the fundamental predictions of RFOT, in particular the 
possibility to force a small enough system into an ideal glass state, and present several new ones, 
concerning aging properties or non-linear rheology. Finally, we compare RFOT to other recent 
theories, including elastic models, Frustration Limited Domains or Kinetically Constrained models. 



Details that could throw doubt on your interpretation must be given, if you know them. You must do the best you 
can - if you know anything at all wrong, or possibly wrong - to explain it. If you make a theory, for example, and 
advertise it, or put it out, then you must also put down all the facts that disagree with it, as well as those that agree 
with it. 

R. P. Feynman 
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I. INTRODUCTION 

When we have to argue why one should be interested in the theory of glasses, it is customary to quote Phil 
Anderson who wrote that the deepest and most interesting unsolved problem in solid state theory is probably the 
nature of glass and the glass transition. He was of course himself deeply involved in this endeavour and as always 
made several outstanding contributions. He was convinced early on that spin-glasses would be a kind of appetiser, 
the theory of which would open up new avenues to understand glasses. As we all know, the theory of spin-glasses 
turned up to be exquisitely complex, offering one of the most beautiful surprise in theoretical physics in the last 30 
years. Very soon after the solution of mean-field spin-glasses was worked out by Parisi, Derrida and others 0, [![, 
Kirkpatrick, Thirumalai and Wolynes, in a remarkable series of papers, proposed the spin-glass inspired 'Random 
First Order Transition' (RFOT) theory of the glass transition [4j, la, la, ll| ■ It took ten more years for the theoretical 
physics community to understand the depth and scope of RFOT, and to start actively criticising, reformulating and 
expanding the original papers and the more recent contributions of Wolynes and associates [1, [l(| [HI • Whereas 
many aspects of the theory are quite compelling, others still appear wobbly. It is fair to say that we are still far from 
a consensus on whether the basic tenets of the theory are correct or not, and whether Anderson's suggestion is the 
right track to follow. 

The aim of this chapter is to summarise the basic arguments and the intuition bolstering the RFOT picture for 
glasses (for recent reviews with a significant overlap, see [HI, O 0] ) ■ We want to review the pros and cons that 
support or undermine the theory, and the directions, both theoretical and experimental, where progress is needed to 
ascertain the status of RFOT. We have deliberately written this paper using a narrative style, insisting on ideas and 
concepts and leaving technical details to original publications or some appendices when we felt they were particularly 
important. We have also chosen to present several conjectures or half-baked arguments that may well be wrong, with 
the hope that these will stimulate more research and discussions. Although we are ourselves somewhat biased in 
favour of RFOT, we are in fact motivated by a strong desire to understand the physical mechanisms underlying the 
glass transition, and hope that the following material will help making headway. 

A. Glasses: back to basics 

So what is so special about glasses, that requires tools unavailable to classical statistical mechanics and solid state 
theory? Coming back to the basic phenomenology, a glass is an amorphous solid: below the glass temperature T g , 
the static shear modulus Go is for all practical purposes non zero (like for a crystal) but there is no apparent long 
range order: glasses are liquids that cannot flow. This is puzzling, because the shear modulus can in principle be 
computed as a purely thermodynamical average of microscopic quantities, since it measures the change of free energy 
to an infinitesimal change of the shape of the container. This is necessarily zero for an ergodic amorphous state 
(liquid) [15lll6j. As Anderson puts it [171] ■ acquiring rigidity is not a minor fact: We are so accustomed to this rigidity 
property that we don't accept its almost miraculous nature, that is an "emergent property" not contained in the simple 
laws of physics, although it is a conseguence of them. 

In the case of a crystal, we understand what is going on: the liquid undergoes a first order phase transition towards 
a symmetry broken state characterised by a periodic arrangement of the particles. The fact that symmetry is broken 
means that ergodicity is broken as well: disordered (liquid) configurations are no longer accessible to the system. 
Glasses, on the other hand, are in an ergodicity broken state but with - apparently - the same symmetry as the 
liquid. Below T g , glasses remain stuck around a mechanically stable, but amorphous configuration that can sustain, 
for a very long time, an external shear without flowing. When microscopic energy barriers A are high, as in so-called 
strong glasses like Si02 , this is simply due to the fact that elementary moves themselves become extremely slow, and 
the viscosity r\ follows an Arrhenius law, r\ = ?yoexp(A/T) which eventually exceeds 10 13 Poise at low temperatures, 
and we conventionally call the system a glass. But in fragile molecular liquids, the energy barrier itself grows as 
temperature is decreased. This suggests that some kind of thermodynamic "amorphous order" propagates (at least 
over medium scales), such that any change of configuration attempting to restore ergodicity requires the collective 
rearrangement of many particles - this is the only way to prevent it from happening very quickly. In other words, the 
energy barrier for the rearrangement of a cluster of particles must grow with the size of that cluster, at least up to a 
certain size. 

From a theoretical point of view, it is convenient to think of an idealisation where the energy barrier and the 
relaxation time are actually infinite, so the system is really in a new thermodynamic phase, and ask oneself how to 
describe this new phase and the transition separating it from the ergodic liquid. After all, this is what happens when 
a system crystallises and one ignores defects, such as vacancies and dislocations. There are several crucial differences, 
though. First, as we alluded to above, there seems to be no simple static order parameter to landmark the transition; 
the structure factor of the glass, for example, shows very little difference, if any, with that of the liquid. Second, there 
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FIG. 1: Excess entropy S xa {T) as a function of temperature for Salol (in kJ/K /mot). The solid line is a fit of the form 
S X3 (T) = A(l — Tk/T). Note (i) the abrupt change of slope at T g , interpreted as the freezing of the configurational entropy 
contribution and (ii) the value of S X s(T g ) sa 3ks per particle. From [181 ]. 



is no latent heat at the glass transition, the glass freezes in whatever entropy was in the liquid at T g , and this is large: 
the excess entropy S xa of the glass over the crystal at the same temperature (in order to remove at least part of the 
vibrational entropy) is substantial, of the order of ~ Iks per molecule [lH - for example, S xs (T g ) w 3ks per particle 
for Salol - see Fig. 1. This means that the number of microscopic configurations in which the glass can get stuck is 
exponentially large. There seems to be a very large degeneracy in the ways molecules can arrange themselves such 
as to form mechanically stable, amorphous patterns around which they vibrate without exploring the other stable 
patterns. As a matter of fact, those thermal vibrations are small, the Lindemann ratio that compares the root-mean 
squa re of the vibration amplitude u to the intermolecular distance a is around 10% at T g for most molecular glasses. 
[214j | Correspondingly, as soon as the glass transition is crossed, the Debye- Waller factor is close to unity. 

This exponential degeneracy of the "phases" towards which the liquid can freeze is the feature that classical theories 
of phase transitions cannot easily handle, and that requires new tools. At the same time, this feature seems to be the 
very essence of glassiness: in order to prevent fast crystallisation, the interaction between molecules must be able to 
generate enough "frustration" to make the energy landscape rough and rocky and trap the system in a configuration 
not very different from an arbitrary initial configuration of the liquid. Before explaining how the solution of some spin 
glass mean-field models has afforded us the analytical tools needed to deal with exponentially degenerate amorphous 
system, we want to delve a little longer into the phenomenology of glass forming liquids, and identify several salient 
observations, common to some degree to all glass formers, and that ought to be explained by any viable theory. 



B. Glasses: more advanced phenomenology 



As recalled above, a distinctive properties of so-called fragile glass formers is the very fast increase of the relaxation 
time (or the viscosity) as temperature is decreased. The slowdown is said to be super- Arrhenius and cannot be 
explained by activation over fixed barriers. It is customary to define the fragility parameter m as: 



d log 



10 'a 



91n(l/T) 



(1) 



where the a relaxation time r a is measured in microscopic time scale (10 -13 sees.). For a purely Arrhenius slowdown, 
the value of m is tuq = 16, by definition of the glass temperature T g . Large values of m — m Q means that the effective 
barrier A(T) that one has to plug into an Arrhenius law strongly varies with temperature. For the most fragile glass 
formers, the effective energy barrier at T g can be 5 times larger than its extrapolated high temperature value. This 
is a very substantial increase pointing towards some collective mechanism. It is clear that many distinct numerical 
functions will do a good job at fitting a five fold increase of an energy barrier over a 30% temperature variation 
(see the discussion in Well known examples of fits are (a) the Vogel-Fulcher form Avf oc TTyp/i/T — Tvf), 

suggesting a divergence of the barriers at a finite temperature Tvf) (b) the Bassler form A# cx T* 2 /T, for which the 
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FIG. 2: Energy barrier A = fcsTlnr a /ro in units of the Goldstein crossover temperature T* , as a function of T* /T, for three 
fragile glass formers (from [23[). One observes (i) the change of regime around T* and (ii) an energy barrier that increases by 
a factor up to five between T* and T g . 



barriers only diverge for T — 21]. T* is the temperature at which the liquid starts to be super- Arrhenius. More 
sophisticated fits, where the temperature dep endent contribution to the effective barrier vanishes above T* , have also 
been proposed in the recent literature [22], [23|, [24| , see Fig. 2. The existence of such a crossover temperature T* , where 
the slowdown mechanism appears to change from weakl y ac tivated to strongly activated has been advocated several 
times in the literature since early insights by Goldstein [251 ]. Glasses lose their rigidity as temperature increases not 
merely because the relaxation time scale decreases, but more radically because at some point local stability is lost. 
In the language of Mode-Coupling Theory [26|, [U [28| , that we will discuss below, local "cages" open when a certain 
critical temperature is reached, above which particles move quasi-freely around one another, and below which thermal 
activation is necessary. We will see below how this crossover naturally appears in the context of RFOT. Fig. 2 in fact 
advocates a cross-over from Arrhenius to super- Arrhenius behaviour at T* , whereas the above MCT interpretation 
suggests that these barriers are actually absent above T*. We will not embark in a detailed discussion of the possible 
origin of this high temperature barrier, but as will be clear below (see in particular section IV C[) . one actually expects 
activation processes to first appear at a higher "onset" temperature To > T*. 

Another striking empirical regularity of glass formers is remarkable correlations between thermodynamics and 
dynamics, for a larg e variety of materials (for introductory material, see [lij [3(J HH H3] ) • These so-called Adam-Gibbs 
correlations [33l.l34| hold both for single liquids as temperature is varied, and cross-sectionally for different liquids. For 
single liquids, it relates the relaxation time to the excess entropy mentioned above, and reads: In r(T) oc [TS^T)] -1 
(see Fig. 3-a). This is the Adam-Gibbs relation. Since the excess entropy appears to extrapolate to zero at the 
so-called Kauzmann temperature Tki the Adam-Gibbs relation is compatible with the Vogel-Fulcher description if 
Tvf ~ Tk- The coincidence between these two temperatures has been documented for many materials [3l| . but has 
also been disputed as an artefact induced by unwarranted fits (see e.g. ]35, 36]). 

Cross-sectionally, one observes substantial correlations between the fragility parameter m and the jump of specific 
heat AC P at T g : m — mo cx AC P (see 37J and Fig. 3-b). More fragile glasses have a stronger variation of their excess 
entropy with temperature and a larger specific heat jump at T g . Note that since super-cooled liquids are formed by 
very different microscopic molecules, it might be reasonable to speak of a AC P per "bead", the beads corresponding 
to the truly mobile elements inside one single molecule Q . This is an approximate way to compare diverse molecules 
in a more meaningful way. With this normalisation the correlation between m and AC P increases substantially [371 ] . 

Although these correlations have been repeatedly reported in a very large body of experimental work, they are 
obvi ously not perfect (see Fig. 3-b) and can always be questioned, especially for super-cooled polymeric melts [11]. 
|215j ] Furthermore, correlation does not mean causality - in particular, it is not obvious that the cause for the super- 
Arrhenius slowdown of glass formers is the behaviour of the excess entropy, a purely thermodynamic quantity. Whereas 
RFOT essentially thrives on such a causal mechanism where thermodynamics drives dynamics, other scenarios have 
been put forth, where thermodynamics is completely irrelevant to understand the dynamics, and any correlation 
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FIG. 3: Left: Illustration of the Adam-Gibbs correlation: plot of log 10 r a /ro as a function of the inverse configurational entropy, 
revealing a reasonably linear region between 10 -7 sec. (corresponding to T = T*) and 10 sec. (T ~ Tg). From Right: 
correlation between the fragility parameter m and the jump of specific heat at T g , AC P , for a collection of 20 molecular glasses 
for which we were able to cross-check the results from different sources. The correlation is significant, albeit far from perfect 
(R 2 = 0.52). The regression gives AC P ~ 0.15m. A better correlation, with much less scatter, is obtained when AC P is counted 
per "bead" [37||. 



between the two is deemed fortuitous. We will come back to this point in sections IVI AIVI Bl 

Finally, as discovered more recently, dynamics is glasses is spatially heterogeneous and temporally intermittent (iol . 
l4lj . A number of experimental, numerical and theoretical papers have established that the slowdown of supercooled 
liquids as temperature is reduced is accompanied by the growth of a purely dynamical correlation length with no 
static counterpart. This length can be for example defined by measuring how a local perturbation of the system (for 
example a local density change) influences the dynamics at a distance r from the perturbation, £<j is the decay length 
of this response function (see [42l. |43|). 

The fact that dynamics in glasses is heterogeneous is thought to have observable consequences on macroscopic 
quantities [4(3, |4lj. For example, the stretched exponential nature e - ' of the relaxation is often interpreted in terms 
of a mixture of exponential functions with different local relaxation time, reflecting the coexistence of slow and fast 
regions. Stronger heterogeneities should lead to smaller values of (3. Similarly, violations of the Stokes-Einstein relation 
(SER) between the viscosity and the self-diffusion constant have also been attributed to dynamical heterogeneities. 
Experimentally, the self-diffusion constant at T g can be 10 3 times larger than expected from the value of the viscosity 
if the SER was valid. There seems to be a strong correlation between the violation of SER and the smallness of 
/3, indicating that both phenomena should indeed have a common origin. Dynamical heterogeneities represent a 
relatively new facet of the glass transition problem. It is certainly a very important one that contribured shifting the 
attention of theorists towards the real space, microscopic origin of glassy dynamics. Their existence demonstrate that 
glass formation cannot be thought of as a purely local process, due to the increase of a local energy barrier, or to the 
decrease of a local free volume, etc... Although dynamic heterogeneity may not be the cause of the slowing down of 
the dynamics, any theory of the glass transition should now account for the rich non-local space-time properties of 
the dynamics. 

Many other interesting phenomena and empirical regularities of glass forming liquids would be worth reviewing, 
but are beyond the scope of the present chapter. Some will however be mentioned below, in relation with theoretical 
predictions. 
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II. GLASS THEORY: WHERE SHOULD ONE START? 



A. General remarks and the Random Energy Model 



As we mentioned above, glasses below T g are for all purposes in a genuine thermodynamic state, characterised by 
a certain specific heat, a shear modulus, a Debye- Waller factor, etc. that one wishes to predict from first principle. 
The fact that after a very long time the system actually relaxes and flow, and is in fact technically speaking a liquid, 
is irrelevant except if one is prepared to wait a million years. After all, even a bona fide crystal subject to shear will 
eventually flow [39(, and diamond eventually turns into graphite. Keeping in mind that the glass transition we are 
talking about is actually waiting time dependent, a natural statistical mechanics approach is to find an appropriate 
mean-field theory that captures the physics qualitatively, even if quantitatively wrong - for example, predicting a true 
transition while there is in fact only a sharp crossover, or predicting spurious metastability wiped out by nucleation 
effects, or inaccurate critical exponents renormalised by strong fluctuations. This is a very usual situation, and we 
know that hard work is usually required to understand and account for non-mean-field effects. The classical example 
is the Curie- Weiss theory of magnets that required more than sixty years and the Wilson-Fisher renormalisation group 
to account for real, three dimensional magnets. 

In the case of glasses, we need a theory with a transition that has a mixed first/second order character. Empirically, 
there is indeed no latent heat nor volume jump at the glass transition, whereas the amplitude of frozen-in density 
fluctuations (again on a long but finite time scale) jumps discontinuously from zero in the liquid phase to a large 
value. As we noted above, as soon as the glass transition is crossed, the Debye- Waller factor is close to unity; in this 
sense the transition is strongly first order. As intuitively expected, it seems impossible to maintain a classical system 
around an amorphous configuration if thermal fluctuations are too strong. 

As anticipated by Anderson, progress indeed came from spin-glasses, albeit in a rather unexpected way. The 
natural mean-field model for spin-glasses proposed by Sherrington and Kirkpatrick (SK) was solved by Parisi, who 
invented along the way a framework ('replica symmetry breaking', RSB) to deal with systems with a large number of 
quasi-degenerate states @j. Whether or not the mean- field limit of spin-glasses provides a reasonable starting point 
to understand real spin-glasses is an open issue which we do not want to even touch upon. What is interesting, in any 
case, is that the SK model is characterised by a genuine transition towards a state for which "long range amorphous 
order" has a precise meaning. Spins freeze in a random configuration but the correlation length of the fluctuations 
are infinite. As one approaches the transition, the correlation length diverges, and some observab le q uantities, such 
as the non- linear (third order) susceptibility, are directly sensitive to this growing amorphous order [44l[45|. However, 
the SK transition is second order, and the spin-glass order parameter (i.e. the frozen-in magnetisation fluctuations) 
is infinitesimal just below the transition, in stark contrast with what happens in molecular glasses. 

Because the mathematics of the Parisi solution is exquisitely baffling, Derrida [|[ sought for a simpler model of 
spin-glasses, and came up with the Random Energy Model (REM), which is at first sight too simple for anything 
interesting to happen. In that model, the energies of microscopic configurations are independent and identically 
distributed (IID) random variables. Still, the thermodynamics is non trivial. There is a second order phase transition 
between a low temperature phase where the extensive part of the entropy is zero and the system is trapped in a 
handful of low-lying energy states, and a high temperature phase where a very large number of configurations are 
relevant. Physically, this comes from the fact that the entropy per spin a as a function of the energy per spin e 
vanishes linearly at the edge of low-lying energy states: 



This assumption is enough to get a transition temperature Tk, since from (dcr/de) = l/T one gets for the entropy as 
a function of temperature S(T): 



cr(e) w l3 K {e- e min ) - — (e - e min ) 2 + . . . 



e > e 



B > 0. 



(2) 




with Tk 



(3) 



and £(T < Tjc) = 0. Clearly, the specific heat has a jump at Tk, given by: 



(4) 



The above shape of S(T) for T > Tk appears to be a good fit of the excess entropy of a number of materials, see Fig. 
1 and [H|]. This, together with the Adam-Gibbs relation lnr cx (TE(T)) -1 , precisely leads to the Vogel-Fulcher law 
for the relaxation time. Note that in the analogy between the REM and real glasses, each REM configuration is an 
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amorphous metastable state in which the liquid can be trapped in. Energies of configurations are of course not IID 
in a real glass but one can argue, as we shall do in the following, that they behave as if they were so. 

The important step forward was to realize that instead of being just an abstract toy model, the scenario laid 
out by the REM is exactly realized for a large family of disordered mean-field models with well defined degrees of 
freedom, which allows one to make bona fide thermodynamical and dynamical first principle calculations [3j, [J, [5|, |46| ] . 
A particularly important model within that family is the 3-spin Hamiltonian: 

H = V3^ ^2 e l jkS l S J S kl S = ±l (5) 

i<j<k<N 

where the sum runs over all triplets of spins and Cijk are IID random variables with zero mean and unit variance. In 
order to compute the thermodynamical properties of this model at low temperatures, one has to use Parisi's replica 
symmetry breaking theory, but in this case the technicalities arc different (and in fact simpler) than for the SK model 
- these models are called '1-step RSB' (1-RSB). One finds a REM type transition with a mixed first/second order 
character. The entropy and the specific heat of the low temperature phase are now non zero but the specific heat jumps 
at the transition temperature Tk- The amplitude of the frozen- in magnetisation fluctuations jumps discontinuously 
from q = for T > Tk to q* ~ 1/2 for T — T^. These similarities with the phenomenology of glasses suggest that 
TRSB models may indeed grasp some essential aspects of the glass formers. This observation motivated a series of 
papers by Kirkpatrick, Thirumalai and Wolynes (KTW), who established deeper still connections [4|, Ifa, [g, [D] , as we 
detail below. 

The only certainty, at this stage, is that there exists well defined, but somewhat exotic mean-field models that 
exhibit a phase transition with a new phenomenology. That these fully connected disordered spin models have 
anything sensible to say about the viscosity of supercooled liquids is obviously highly questionable, but it was the 
hunch of KTW that they do. 



B. Physical properties of 1-RSB models 

1 . Thermodynamics 

1-RSB models are disordered models solved by the so-called one-step replica symmetry breaking ansatz (see e.g. 
[Ill IS Eil and references therein). Examples are the 3-spin model defined above, or any p-spin model with p > 3 (the 
REM corresponding mathematically to p — > oo 0]). A 2-spin model with 1-RSB properties is the Random Orthogonal 
Model where the interaction matrix Jy between Ising spins is a random matrix with all eigenvalues equal to +1 or 
— 1 [HoL IHH l52| . Other examples of direct physical interest will be given below. 

The physical behaviour of these models is particularly transparent in terms of the Thouless-Anderson-Palmer 
(TAP) 53| approach that allows to scrutinise the free energy landscape of the model. Technically, one computes the 
Legendre transform of the free energy as a function of all the local magnetisations. The minima of this TAP free 
energy correspond to the thermodynamic states of the system, like the two minima of the Curie- Weiss free energy 
represent the two low temperature ferromagnetic states in the Ising model. The configurational entropy, or complexity 
£(T) is defined as N -1 In Af where N is the number of solutions of the TAP equations that contribute to the free 
energy density. In the limit p — > 00 the TAP states of the p-spin model become the REM configurations. 

For a TRSB system, the analysis of the complexity shows that on top of the temperature Tk below which E(T) 
vanishes, there exists a second, higher temperature that we call Td > Tk (the subscript d will become clear below) 
above which E(T) drops discontinuously to zero again. However, the situation is markedly distinct below Tk and 
above Td- As these two temperatures are crossed, the properties of the free-energy landscape change drastically, but in 
two very different ways. At Tk the number of minima is no more exponential in the system size, X(T < Tk) — 0, and 
this leads to a thermodynamic transition of the REM type, as explained above. The only difference with the REM is 
that now the free energy minima are not single configurations; they also involve small vibrations around the minima 
that contribute to the total entropy in this phase. At Td, on the other hand, there is no thermodynamical phase 
transition at all, but rather a "fragmentation" of phase space into disconnected minima. For T > Td there is only one 
minimum of the TAP free energy with / = f P (T), where all local magnetisations are zero, corresponding to a trivial 
paramagnet. For Tk < T < Td, on the other hand, the partition function is built up from the superposition of an 
exponentially large number of disconnected minima, which all have the same free energy / = f*(T). As a consequence 
the free energy of the system, still given by f p (T), is smaller than f*(T) because of the additional entropic gain due 
to the complexity: 



f p (T)=f*(T)-TE(T), 



(6) 
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FIG. 4: Summary of the mean-field 1-RSB scenario for p-spin models: free-energy of the different TAP states as a function 
of temperature, from the lowest (most stable) f m in(T) to the highest (marginally stable) threshold states. We also show the 
free energy of the paramagnet f P (T) that crosses / m in(T) at the static transition Tk- We also show the free energy f*(T) 
of the dominant TAP states in the temperature range T £ [Tfc,Td], and the configurational entropy E(T) that makes up the 
difference between /* and f p . The most stable branch terminates at the onset temperature To. 



However, this result implicitly assumes that all these minima are mutually accessible, so that the corresponding 
subtraction of T£(T) is warranted. Quite surprisingly, one finds that f p (T) is not singular at T^: the thermodynamics 
is completely blind to this fragmentation! We will come back to this crucial point in the following subsection. 

The abrupt vanishing of £(T) above Td does not mean that there are no minima in the high temperature region, 
but rather that these states are no longer numerous enough to compete with the trivial TAP minimum where all local 
magnetisations are zero. More formally, one can introduce a free-energy dependent complexity, a(f,T), that counts 
the nu mber of TA P minima with free-energy density / at temperature T. The partition function of the system then 
reads: HH EH 



Z(T) = J df exp 



(7) 



For large N, one can as usual perform a saddle-point estimate of the integral, that fixes the dominant value of /, 
f*(T), such that: 



T da(f,T) 



df 



= 1 (8) 

f=f'(T) 



The temperature dependent complexity is in fact defined by: £(T) = o~(f*(T),T). From this point of view, the 
fact that f p (T) is not singular at Td is actually far from obvious (see Fig. [3}. It is one of the most unexpected 
result emerging from the analytical solution of 1-RSB models. It suggests that at Td the paramagnetic (or liquid) 
state fractures in an exponential number of states and that this transition is only a dynamic phenomenon with no 
consequences on the thermodynamics. This physical scenario is key to the development of the RFOT theory of glasses. 
When T > Td, there are still non-trivial TAP states but their contribution is sub-dominant with respect to the 'trivial' 
paramagnetic state, which exists in that region. The probability to sit in one of these TAP states is exponentially 
small in N ior T > Td. Technically, this means that a(f*(T),T) becomes too small to compensate the difference 
between /* and f p . One has reach a third temperature To > Td for the last TAP state to disappear, in the sense that 
a(f, T > Xb) = 0. A sketch of the rather intricate situation is provided in Fig. [U The temperature T turns out to be 
important for the understanding of glasses, because the existence of locally stable states can slow down the dynamics 
of the system. It is often called the 'onset' temperature, where activation effects first appear. 

Another way to describe the above fragmentation scenario is to compute the free energy cost V(q) to keep an 
arbitrary configuration of the spins {Sf}, drawn from the canonical probability distribution at temperature T, at a 
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certain overlap q = -A- ^ S^S? from another equilibrium configuration {Sf } The form of V(q) for T > Td is 

such that the only minimum of V(q) is at q = 0, for which V{q = 0) = 0: spontaneously, two copies of the same 
system wander in phase space independently and the most probable is to find them in completely different states. 
When Tk < T < Td, a second minimum appears for a value of q* > 0, but the absolute minimum is still at q = 0, 
and V(g*) > is precisely equal to the complexity S(T). This is perfectly in line with the above interpretation: the 
secondary minimum corresponds to both systems being in the same locally stable TAP state. However, the associated 
loss of configurational entropy makes this situation metastable, and the two copies will end up in equilibrium in two 
different states with zero mutual overlap. Finally, when T = Tk, V(?*) = and there is a finite probability to find 
two independent copies of the system in the same state. 



2. Energy landscape 

It is useful to rephrase, and make more precise, the above discussion of the free-energy of the TAP states in terms of 
the properties of the energy landscape. The topological properties of the Hamiltonian of the system as a function of the 
microscopic configurations are obviously temperature independent. Particularly important configurations are those 
corresponding to local minima. The TAP states can be seen as the ensemble of thermally populated configurations 
around a local minimum. The spread of the statistical weight around the minimum gives rise to magnetisation 
fluct uations and to a vibrational entropy contribution s v , which depend on the local curvature of the Hamiltonian. 
[218j | Each TAP state is therefore characterised by a 'bottom-of-the-well' energy eo (per spin) and an entropy s v , giving 
rise to a free energy / = eo — Ts v . In the low temperature harmonic approximation, the fr ee-en ergy per particle can 
be expressed in terms of the density of eigenvalues p(X) of the Hessian stability matrix as: |219l ] 

/(T) = e + | J dXp(X) In A - | 1h(2ttT). (9) 

It turns out that for p-spin models, the density of eigenvalues is a shifted Wigner semi-circle, 



/w ^ v(A — A min )(A max — A) 
P A) = 7T V2 , (10 

7T (A max - A,, ' '- 



v min j 



with a lower edge of the spectrum A m i n > that only depends on the bottom of the well energy eo [55| • The remarkable 
feature is that A m i n decreases monotonously when eo increases, i.e. lower energy states are also more stable in the 
sense that the average curvature of the Hessian is larger, and thermal fluctuations weaker. There exists a particular 
'threshold' value eth of eo at which A m ; n is exactly zero, i.e. states with the threshold energy are marginally stable. 
There is an exponentially small number of minima with energy above eth , stationary points of the Hamiltonian are 
mostly unstable saddles with a negative A m i n ; |A m i n | in fact increases linearly with eo — eth close to the threshold. 

Interestingly, one finds that the bottom-of-the-well energy of the TAP states that are dominant right at Td is 
precisely eth (see Fig. 01 and therefore these TAP states are also marginally stable and have soft modes. This 
marginality has important physical consequences, in particular because it induces the divergence of a correlation 
length close to Td (see Eq. (|25p below and the discussion there.) Note also that energy minima with eo = e t h are 
exponentially more numerous than any other minima. This means that the overwhelming majority of minima of 
1-RSB systems are only marginally stable! 

On the other hand, the bottom-of-the-well energy of the TAP states that are dominant at Tk is clearly the ground 
state energy e m ; n , and these states are also the most stable ones, with the largest value of A m i n . 

As the temperature increases, the harmonic approximation breaks down and at some temperature TAP states "melt" 
and disappear (see Fig. 2]). This melting temperature depends on eo, and (as intuitively expected) increases as eo goes 
down, since the stability of the corresponding states increases. This is strongly reminiscent of the recent experiments 
on ultra-stable glasses by Ediger [56|. Furthermore, as we shall discuss later, TAP states that are irrelevant for the 
thermodynamics within mean-field become important in finite dimensions, especially for the dynamics that becomes 
activated above Td- From this perspective it seems natural to associate the onset temperature To mentioned above 
to the melting of the most stable TAP states with eo = e m i n [57j ■ 

The interest of exactly solvable models in general is that they allow to characterise the physical properties of a 
system in a solid, rigorous way. These solvable models sometimes unveil effect that cannot be easily guessed using hand 
waving and plausibility arguments. In the case of p-spin models, analytic calculations indeed provide an extraordinary 
rich information on the structure of the energy landscape and the thermodynamical properties, which we tried to 
summarise in the above two sections. The most striking and unexpected result is the fragmentation transition at 
Td, where the nature of paramagnetic phase changes without any thermodynamical signature. This phenomenon is 
bound to have some impact on the dynamics of these TRSB models, to which we turn now. 
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3. Dynamics 



Below Td, stable local minima are numerous enough to dominate the thermodynamics. These local minima have 
typically zero mutual overlap, which means that the local amorphous order is very different in the different TAP 
states. But since we are in mean-field, the energy barrier between these zero overlap states involve all spins and 
must diverge in the thermodynamic limit N — > oo. Dynamically, the system must thus be trapped in a single TAP 
minimum when T < X^; the thermodynamical average is only realized if the large time limit £ — > oo is taken before the 
limit N — > oo. The temperature Td should therefore separate an ergodic high temperature region from a glassy phase 
disguised in a 'non ergodic paramagnet' below Td- These intuitive expectations are confirmed by the exact treatment 
of the dynamics of these 1-RSB models, assuming Langevin dynamics and continuous spins. The dynamical equations 
can be cast in terms of two observables: the spin-spin correlation function, defined as: 

i 

and the corresponding response function R(t) = d(Si(t)) /dhi(Q), where /ij(0) is the local magnetic field at site i and 
time 0. In equilibrium, the two functions are related by the Fluctuation-Dissipation Theorem: TR(t) — —dtC(t). In 
the large N limit, and for the 3-spin model, C(i) is found to obey a closed non- linear integro-differential equation 
that reads (see jllilloll): 

d t C(t) + Tv l C(t) = J d S C 2 (s)d t - s C(t - s). (12) 

Remarkably, this equation is identical to the so-called schematic Mode-Coupling approximation, developed since the 
mid-eighties by Gotze and others as a self-consistent re-summation scheme for the density correlation function of 
strongly interacting liquids [2(| H?], [H, [59J. The solution of the above Mode-Coupling equation is now very well 
understood mathematically. At high temperature, interaction is unimportant and C(t) relaxes exponentially. As the 
temperature is reduced towards Td = y / 3/8Jo, the relaxation becomes slower and non-exponential. There is a first 
rapid decay from C(t = 0) = 1 towards a plateau value q* , then a slow evolution around it (the (3 regime) and finally 
a very slow decay towards zero (the a relaxation). At a given temperature T > Td, the dynamics in the /3-regime is 
described by power-laws. The approach to the plateau can be written as C(t) ~ q* + ct~ a , and the later departure 
from it as C(t) ~ q* — c't b , where c, c' are constants, and a and b are exponents, found to be a rj 0.395 and 6=1 
for the 3-spin model. In the a-regime close to Td, C(t) verifies a scaling law, called Time- Temperature superposition 
(TTS) in the structural glass literature: 

C{i) ~ g(t/T a ) , (13) 
where g is a certain scaling function and r Q is the a relaxation time that diverges at the transition T = Td as: 

r a( x(T-T d )-^ 7 =_L + i. (14) 

la lb 

When T < Td, the asymptotic limit of C{t) is non longer zero, which means that ergodicity is broken and the 
assumptions made to derive Eq. (|12p are not valid. For the 3-spin model, one can still write down the dynamical 
equations in terms of two-time (aging) quantities, C(t, t w ) and R(t, t w ) - see [H,[63| for details. If the initial condition 
is taken instead with the Boltzmann measure at a temperature T < Td (6l| one finds a rapid decay of C(t) towards a 
non zero value. This corresponds to the equilibrium relaxation inside one of the TAP states with free energy f*(T). 
No paramagnetic or liquid state is found below Td within this approach, in agreement with the fact that it does not 
exist anymore since it is fractured into an exponential number of TAP states. 

These findings are perfectly in line with the energy landscape picture brushed in the previous paragraph (see 
e.g. [48|, 162J, |63|): for T > Td, the system only explores with high probability unstable stationary points of the 
Hamiltonian, and relaxation can proceed without any large barrier crossing: the correlation function decays fast and 
quasi-exponentially. As T approaches T^, both the number of unstable directions and the negative curvature of these 
unstable directions go to zero. One therefore expects a fast relaxation to a plateau, governed by the stable directions, 
followed by a slow relaxation driven by the unstable directions that allow the system to relax towards zero (see [64| 
for a simple picture, and [65[ for more insights). Since the curvature of these unstable directions goes to zero, the 
relaxation time diverges continuously at Td- Below Td, all directions are stable and the system is trapped within a 
TAP state. The only relaxation mechanism is (collective) activation, but for mean-field models with N — > oo this 
relaxation channel is also forbidden, and Td is a genuine ergodicity breaking solution, precisely of the same nature as 



13 



found within the Mode-Coupling Theory (MCT) of glasses. Although it is very clear that non mean-field activated 
events will smear out this dynamical transition, the precise theoretical description of this mechanism is well beyond our 
current abilities. As we will detail below, the RFOT framework only brings partial answers to this crucial question. In 
any case, the similarity between this smeared transition at Td and the dynamical crossover conjectured by Goldstein 
around T* is very striking. 

It should be noted that the above discussion is not at all limited to the 3-spin model. All models that exhibit a 
1-RSB transition and for which the dynamics can be studied analytically share the very same features. Although 
the corresponding integro-differential equation for C(t) can take more complicated forms, the two-step relaxation 
described above, with power-law dynamics around the plateau and a diverging a relaxation time. Only the value of 
the plateau q* and of the exponents a and b are model dependent (although the relation between a and b is not). In 
fact, the MCT phenomenology is even more general since all the results can be derived based on general assumptions 
about the nature of the dynamical arrest without relying on any particular model [6t| . 

C. But are we on the right track? 

Let us summarise where we are: we have a mean-field model of spin glasses with random multi-spin interactions. 
This model has a static transition temperature Tk with a vanishing excess entropy a la Kauzmann and a dynamical 
MCT transition at a higher temperature T* ~ Td below which thermal activation becomes dominant, a la Goldstein. 

The analogy may however look far-fetched or even dubious to many. According to Jim Langer, for example, it is 
hard to believe that long range interacting spin model with weird interactions can teach us anything about molecules 
that stop jittering around [671 ]. There are indeed different issues that need clarification. First, in what sense a model 
with spin degrees of freedom on a lattice can be used to describe positional degrees of freedom of molecules? Second, 
spin- glass models assume from the start the existence of quenched disorder and encode it as random interactions, 
whereas molecules in a glass interact via rather simple deterministic potentials. Disorder in the latter case is self- 
induced, in the sense that the system freezes in an amorphous state, so that the potential seen by a given molecule 
can indeed be considered as random. But the freezing phenomenon is precisely what we want to model! The initial 
assumption of quenched randomness put by hand may not be warranted. Third, non mean-field, finite dimensional 
effects must play an important role. As we just discussed, the MCT transition must be smeared by activation events 
and is at best only a crossover in finite dimensions. Does the physics of the 1-RSB transition survive in finite dimension 
or is it totally wiped out by fluctuation effects? Even more fundamentally, we have to address Langer's strongest claim 
[6?J : I do not think these mean-field models can be used to predict the divergence of the viscosity. On the contrary, I 
think the mechanisms that produce molecular rearrangements in glasses must be localised and that long-range models 
inevitably fail to describe such mechanisms properly. In the next subsections, we try to tackle these three issues. 

1. Mean-field models with spins and disorder? 

The first two points are actually deeply intertwined, and the answer seems to be that the 1-RSB scenario established 
within some specific spin-glass models is in fact far more general. The situation is analogous to the standard Curie- 
Weiss scenario for the ferromagnetic transition, or the Van der Waals theory of the liquid-gas transition, which both 
end up being very generic. There are many complementary ways to show this result. First, although historically 
researchers have started to focus on mean-field disordered spin systems, it is now clear that the mean field theory 
of bona fide schematic models of glass-forming liquids leads naturally to the 1-RSB universality class. An example 
is provided by the so called lattice glass models [68|, [69|, [7(|. These are hard particles lattice models (devoid of 
any quenched disorder) that do reproduce the phenomenology of glasses [7lj and, when solved within the Bethe 
approximation, all display the 1-RSB physics described previously. 

Following a completely independent path, variational or numerical treatments of the Density Functional Theory for 
liquids have been performed in order to study whether amorphous solutions, akin to the TAP states appearing below 
Td, emerge at low temperature or high density [72l [73|. The answer is, again, positive: there are indeed amorphous 
density modulations, which are solutions of the variational problem below a certain temperature. The variational 
parameter is the width w of the localised density peaks around an amorphous packing; the liquid corresponds to this 
width going to infinity. The free-energy as a function of the inverse width F(w~ 1 ) looks exactly as the potential V(q) 
described in section Hi B 1[ with the same evolution as a function of temperature. As a matter of fact, q and w _1 play 
the very same role. 

Another independent route that has been followed to show the connection between 1-RSB physics and glass-formers 
is a formulation a la Landau of the problem. It has been shown that natural approximation schemes for supercooled 
liquid dynamics lead to MCT equations, which mathematically describe the progressive fragmentation of phase space 
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into local minima. As recently argued in 66], MCT can be indeed formulated as a Landau theory where the "order 
parameter" is the difference between the time dependent correlation function and its plateau value. 

There is another, perhaps deeper interpretation of the success of 1-RSB theories, based on the idea of generic 
properties of large dimensional energy landscapes. The energy landscape is a 'height' function that gives the energy 
as a function of the N degrees of freedom of the system under study. It is made of valleys and peaks and saddles, 
and it is natural to ask about the statistical topography of such a landscape (80j |. such as the number of minima or 
maxima around a certain level, or the index distribution of the saddles (the index is the number of unstable directions 
of the Hessian matrix), etc. It turns out that for Gaussian random landscapes in large dimensions, the replica theory 
becomes exact and allows one to classify all landscapes into two distinct categories. One, corresponding to short- 
range correlated landscapes, is precisely described by the 1-RSB phenomenology detailed above. In particular, one 
observes a de- mixing phenomenon where (nearly) all saddle points are minima below a certain level e t h, whereas 
(nearly) all saddles above e t h are unstable [il. l55l l79l]. The second class, corresponding to long-range correlated 
landscapes, requires a 'full' RSB solution analogous to re gular spin-glasses with a continuous transition and describes 
a hierarchical landscape with valleys within valleys, etc. [2201 ] If the energy landscape is not a random function, but 
a generic deterministic function of all its arguments (i.e. the position of all particles in the case of glasses), it may 
be reasonable to imagine that the above classification still holds, exactly as predictions from Random Matrix Theory 
accurately describe the spectral properties of generic deterministic matrices. Up to now, all reasonable approximations 
applied to models that are 'complex enough' to be glassy have indeed led to a 1-RSB transition. An interesting explicit 
case is the frustrated Coulomb model studied by Schmalian & Wolynes [H[ , see also 0, HH [86[ for other similar 
inve stig ations. In fact, the energy landscape of many hard satisfaction problems is exactly of the 1-RSB nature: see 
[zljUJ for more on this aspect. 

Finally, let us illustrate the idea of 'self-induced' disorder on the example of a ferromagnetic 3-spin model on a tree, 
where the J^/. appearing in Eq.([5]) are actually all identical for triplets of spins that are connected. It turns out that 
a numerical simulation of the dynamics of the model starting from a random configuration leads to very strong glassy 
effects, with the same phenomenology as the disordered 3-spin model [8a ]. This is also confirmed by a theoretical 
calculation, where a 1-RSB transition is again found. Intuitively, this comes from that fact that JSiSjSk can be seen 
as a two-spin interaction model between j and k with a configuration dependent, random interaction J^iSi. It 
might be that this idea is much more general (see e.g. [5(1 [H, US HH ) , and that it is possible to coarse-grain models 
of interacting particles in such a way to generate these higher-order, effectively disordered interactions. 

2. Replicas for molecules? 

As it should be clear from the discussion above, if we remove all these historical scaffoldings, it is now indisputable 
that there are very deep connections between the 1-RSB universality class and the glass transition, although one should 
of course worry about how fluctuations effects not contained in mean field theory could change or even completely 
wipe out the mean-field physics. One of the most relevant theoretical breakthroughs allowing one to get quantitative 
predictions using these ideas is the successful transplant of the replica formalism to describe systems without disorder 
(see the contribution of Mezard and Parisi in this book 14]). A very important aspect of the theory, is once again the 
emergence of a 1-RSB solution to describe systems with a large number of metastable states: we show this in detail in 
Appendix A. Some general results on the glass transition, free of any approximation, can in fact be obtained within 
this approach (see Appendix A). The development of replica theory opened the path to a series of new analytical 
approximations to calculate the phase diagram and thermodynamic properties of realistic models such as hard- and 
soft-sphere systems, binary Lennard-Jones mixtures, etc. in the glassy region. Several predictions based on the 1-RSB 
formalism turn out to be in remarkable quantitative agreement with numerical simulations, including the existence 
of two transition temperatures Td and Tk, the temperature dependence of the configurational entropy, and a jump of 
specific heat at Tjc (see [13] )• The quality of the predictions are particularly impressive for the case of hard-spheres, 
and the replica formalism has lead to new insights on that problem, in particular regar ding the value of the density 
at Random Close Packing and the existence of a unique Jamming density - we refer to [lj, [75|, [76| for details. More 
recently, Yoshino and Mezard [77j have shown how to compute within the same formalism the value of the static 
shear modulus Go = G(lu — > 0), that we have argued in the introduction to be the most distinctive feature of the 
glass state. As expec ted on general grounds, one finds the thermodynamical shear modulus to be zero above Tk and 
non zero below Tk- [22 1| 

The reasons why the replica formalism is tailored to capture the phenomenology of glasses are discussed in the 
Appendix A, and in full detail in [l4|. In a nutshell, this comes from the fact that there is no a priori way of 
introducing an external field that selects the amorphous low energy configurations relevant in the glass phase. The 
trick [7^] is to introduce m — 1 additional 'clones' of the same system with a small attractive interaction between 
them, that plays the role of a self-adaptive field guiding the system into the 'right' amorphous configuration (see 



15 



Appendix A and [lij for more precise technical statements) . If the m clones have a non zero probability to end up in 
the same configuration even when the attraction goes to zero, the system is a glass; if the clones all choose different 
configurations, the system is still a liquid. Quite remarkably, the computation of the 'cloned' partition function allows 
one to access the configurational entropy of the system, see [3, [78| and Appendix A. 

Several quantitative predictions can therefore be obtained both for the dynamics (within MCT) and for the statics 
(using replicas). Reviewing them all is beyond the scope of this chapter (see PJj]). However, it is important to 
underline that the possibility of obtaining several quantitative predictions is certainly a very strong advantage of 
RFOT when compared to other theories that often only provide qualitative insights. 

3. Finite dimensional effects? 

How much of the 1-RSB phenomenology survive in non mean-field, finite dimensional situations? Some authors, 
like Mike Moore, have argued that 1-RSB models are inherently unstable in finite dimensions, and fall in another 
universality class, that of p = 2 spin-glasses in a field, for which the transition in continuous or in fact a mere crossover 
[92l [93j . If this were the case, then all the interesting results quoted above would obviously be useless to account 
for the properties of real glasses. This worry stemmed in part from the fact that all spin (or Potts) models with 
a 1-RSB transition in mean field behave totally differently when simulated on a finite dimensional lattice [li], [95[ . 
All MCT/1-RSB features (such as a two-step relaxation function) seem to disappear and leave way to a continuous, 
spin-glass like behaviour. 

Eastwood & Wolynes [l(| have claimed that this is because these models are not "hard" enough (in the sense of 
the height of the non ergodic parameter q* or the Lindemann ratio) to sustain a discontinuous transition in finite 
dimensions. It would be interesting to have a precise formulation of their heuristic argument, and decide from first 
principle whether or not a given model exhibits the 1-RSB phenomenology in finite dimensions. Fortunately, some 
models that are described by 1-RSB in mean field have recently been found to maintain the same behaviour in finite 
dimensions behaviour, at least above the dynamical/MCT temperature T4 - for example the Random Orthogonal 
Model [9^ or several lattice glass models [lH, 197} ■ These findings are important because they shows that there does 
not seem to be any general principle that rules out the relevance of the 1-RSB scenario in finite dimensions, at least 
over experimentally relevant time scales. 

However, one should indeed expect drastic changes going from mean-field models to finite dimensional systems. 
We have already mentioned that activation effects must smear out the dynamical/MCT transition that takes place 
at Trf. As we will argue in the next paragraph, the multiplicity of states with infinite lifetime found in mean-field 
theory cannot persist in finite dimension. Understanding the mechanisms leading to relaxation requires therefore a 
correct treatment of the inter-state dynamics. Make no mistake, Langer is right: these mechanisms, which produce 
molecular rearrangements in glasses, are indeed localised in space. A real-space understanding of finite dimensional 
1-RSB physics in mandatory. But this is precisely the motivation of RFOT! We now turn to this crucial development 
of the theory. 

III. FINITE DIMENSIONS: DROPLETS, CAVITIES & RFOT 

A. The original nucleation argument 

As we have seen above, the dominant TAP minima at temperature T €]Tif,Ty have a free-energy /* per particle 
that is strictly larger than the paramagnetic free-energy f p at the same temperature. Such a situation also arises 
when a ferromagnet in a non zero magnetic field h > is studied in mean-field. When h is not too large, one finds 
that the free-energy per spin has two minima, /+ and /_, corresponding to positive and negative magnetization per 
spin, with /_ > / + . The negative magnetization state is metastable, but the favoured up state cannot nucleate in 
the mean-field limit. In any finite dimension d, the well-known argument giving the critical nucleation radius is to 
compare the free-energy gain of a droplet of up spin in a sea of down spins. This balance reads: 

A(R) = (f + - f_)n d R d + rs d R d -\ (15) 

where Sd, &>d = Sd/ d are respectively the surface and the volume of a d-dimensional sphere with unit radius, and T the 
surface tension between the up and down states. The maximum of A(R) is reached when R = R* = (d—l)F/ (/_ — /+). 
Droplets with size R < R* shrink back to zero, while droplets with R > R* are unstable and grow indefinitely (see 
e.g. [98| for a good introduction). The time to homogeneously nucleate a marginally unstable droplet of size R* is 
t(R*) oc exp(A(i?*)/T), which can be very long when (/_ — / + ) ~ h — ► 0. In the mean-field limit d — > 00, R*(d) 
diverges and one indeed finds strict metastability in this case. 
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In their early paper [g], KTW argued that a similar mechanism is at play for 1-RSB models in finite dimension, 
with /_ — /_)_ = /* — f p = TE(T). In this case, they insisted, the mechanism driving nucleation is of entropic origin, 
which is the difference between one single TAP state and the liquid state. The only difference with nucleation would 
come from an anomalous dependence of the interface energy on i?, which KTW postulated to be of the form Tq(T)R 6 , 
where 9 is a certain exponent < d — 1 and Tq(T) a temperature dependent effective surface tension (that contains 
a d dependent prefactor akin to Sd, which diverges when d — ► oo). Following the same train of thoughts, one would 
conclude that a given TAP minimum is unstable against the nucleation of a paramagnetic droplet of size larger than 
a certain £*, given by: 



dVL d TY,{T) 

However, in the mean-field model, the paramagnet is nothing but the superposition of all relevant TAP states, so it 
is unclear in which sense this whole superposition can "nucleate" as a standard droplet. Furthermore, if this were 
the case, it would mean that TAP states can only be relevant in a transient, non equilibrium regime - once the 
paramagnet has nucleated and invaded the whole system, TAP states would disappear, end of story. And if it is 
one of the exponentially many other TAP states that nucleates, why is the configurational entropy a driving force, 
since one particular TAP state has been singled out? To say the least, the original KTW argument left a many 
people flummoxed and unsatisfied for a long time. Still, their interpretation of the mean-field scenario in terms of a 
"mosaic" of local TAP states with a typical size equal to £* is, we believe, correct, and becomes quite compelling of 
one rephrases the argument as follows. 



B. Entropy driven cavity melting 



We have argued that the 1-RSB scenario is expected to hold not only for abstract spin models but also for realistic 
particle models, on which we now focus. In order to interpret this scenario in finite dimensions, we consider the 
following Gedanken experiment (99l | . Suppose we can identify one of the exponentially numerous TAP state relevant 
at a given temperature T, which we call a, and characterise the average position of all the particles in that state. 
We will establish that there exists a length scale above which the assumption that this TAP state has a well defined 
meaning is inconsistent. In order to do this, we freeze the motion of all particles outside a spherical cavity of radius 
R and focus on the thermodynamics of the particles inside the sphere, S(R), that are free to move but are subject to 
the boundary conditions imposed by the frozen particles outside the sphere. Because of the 'pinning' field imposed by 
these frozen particles, some configurations inside S(R) are particularly favoured energe tically. When T,(T)R d is much 
larger than unity there are many metastable states accessible to the particles in the cavity. The boundary condition 
imposed by the external particles, frozen in state a, act as a random boundary field for all other metastable states 
except a itself, for which these boundary conditions perfectly match. Any other metastable state 7 has a positive 
mismatch energy, otherwise our assumption that state a is locally stable would be violated. We assume this interface 
energy can be written as T a .^R , where T Q;7 > is distributed around a typical value equal to To and 6 < d — 1. 
For simplicity, we assume in this section that T Q;7 does not fluctuate much and is equal to To (see next section for 
a discussion of the effect of fluctuations). We first imagine that we wait long enough so that the cavity embedded in 
state a is fully equilibrated. The partition function Z a can then be decomposed into two contributions: 



Z a {R,T) = exp[-n d R d ^} + ex P 
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(17) 



where / 7 is the excess free energy per unit volume of state 7. (Here and below, all lengths are in units of the 
inter-particle distance a). When T — > Tk, the relevant R is large, allowing one to make clear-cut statements using 
saddle point arguments. We focus on a typical state a at that temperature, i.e. a state with the free energy /* that 
dominates the integral over / above, such that Tda/df = 1. The partition function of the cavity immersed in the a 
state becomes independent of a and reads: 



Z(R,T) 



exp[-ft d i? d ^] 



1 + exp[fi dJ R d £(T) - 



T Q R e 
T 



(18) 



The above expression is central to our argument. When R is smaller than the length £* defined by Eq. ([TBI) but 
still large, the second term is exponentially small even if an exponentially large number of terms contribute. The 
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mismatch energy dominates and the state a favoured by the boundary conditions is the most probable state, even if 
the particles inside the cavity are free to move. In this sense, the cavity is in a glass phase, where only one (or a few) 
amorphous configurations, selected by the boundary conditions, are relevant. When R > £*, on the other hand, the 
second term becomes overwhelming. There are so many other states to explore that it becomes very improbable to 
observe the a state. There leads to "entropic melting" of the cavity. 

In order to derive Eq. (|18|) . we have assumed that the free energy displays many minima, which is only true in 
mean-field. Although in finite dimension the free energy has to be convex, the above procedure is correct: as for usual 
nucleation theory, one starts from an approximate (mean-field) expression for the free energy with many minima, then 
account for the fluctuations that make these states unstable, sum over these states and eventually obtain a convex 
'true' free energy. However, it is instructive to reformulate the previous results directly in terms of the true TAP 
free energy F. Since the cavity is a finite dimensional and finite size system, the free energy must be convex with 
a single minimum. The existence of the cross-over length £* means in this context that for all lengths R < £* there 
is a set of very different boundary conditions indexed by a such that F has a single minimum which corresponds to 
a density profile which is very close to the a one, as obtained from the approximate TAP free energy. However, for 
R > £* , the density profile at the minimum becomes insensitive to the boundary conditions, whatever they are, and 
is a structure-less liquid-like profile. This teaches us that TAP states are, strictly speaking, only defined on scales 
R < £* . On larger length scales the system must "phase separate". In the present case this phase separation is 
not macroscopic as in the liquid-gas problem, [2221 ] but microscopic. This is the essence of the mosaic state that we 
explain in detail below. Note that in practical applications, obtaining a TAP state with the procedure outlined above 
is impossible. A natural approximate procedure, which is commonly used, consists in associating TAP states with 
inherent structures. Although this is quite reasonable, one should not forget that energy and free-energy mi nima are 
not the same concepts at all jlOll Il02j | otherwise paradoxes emerge (like Tk being necessarily equal to zero [103| |). 

C. Properties of the mosaic state 

The argument above indicates that the notion of a TAP state a can only be self-consistent if one restricts to a small 
enough region of space, such that the state is stable against spontaneous fragmentation within its bulk. As we just 
discussed, if one considers true TAP states then there is only a single minimum: the liquid state, independently of the 
boundary conditions above £* and essentially the amorphous state selected by the boundaries below £* . Interestingly, 
cavities smaller than £* behave effectively as if in a true thermodynamic glass phase, even when T > Tk. If on the 
other hand one considers approximate (finite lifetime) TAP states, the TAP states above £* should be obtained as a 
"product" of those obtained on scale £*. As a consequence the free energy of the system should be computed using 
Eq. (JTSJ) above with R ^> £* . This is indeed precisely equal to that of the liquid, where the configurational entropy 
is subtracted from /*. Up to a sub-leading excess contribution from 'domain-walls', one finds: 

fu g (R » £*) « f*{T) - T£(T) + 0{t e - d ), (19) 

in close analogy with the mean-field result, Eq. ^ above. 

The resulting state is called a mosaic liquid. How should one interpret this? Thcrmodynamically, the mosaic 
state is a superposition of all possible TAP states with free-energy /*, as we just discussed. But a reasonable physical 
interpretation is that it is an ever evolving patchwork of "glassites" (by analogy with "crystallites" in a polycrystalline 
phase) of size ~ £*. The bulk of the glassites have a free-energy density /* (up to fluctuations of order £*~ d / 2 , see 
section UlI Fl below) . whereas some excess energy of order £* s ~ d is localised within the grain boundaries between these 
glassites. One expects the local mean-squared particle fluctuations (u 2 (r)), computed over a time much longer than 
To but much shorter than r a , to reveal spatial correlations on a scale £*. Curiously enough, we are not aware of any 
numerical work systematically testing this idea. Appendix B develops and makes more precise these ideas on the 
simpler, but exactly soluble case of the 1-d Ising model. 

Several additional comments are in order: 

• When £* becomes large, the time needed to equilibrate the cavity will also be large, and the above thermo- 
dynamical computation may not be warranted. We will discuss these dynamical aspects in a separate section 
below, but already note that below T g , the system is out of equilibrium and the characteristic size of the mosaic 
is basically stuck at £*(T g ). The free energy of the glass is therefore: 

f glass (T)nf g (T)-TZ(T g ), (20) 

where f g (T) is the free-energy of the most probable states at T gi i.e. f*{T g ) = f g {T g ) (see Fig. [4} There is 
no entropy discontinuity at T g , but a jump of specific heat given by the sum of two contributions: a) since the 
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configurational entropy is frozen, one contribution is AC P = T g dT,/dT\T g , and b) since the local TAP states 
below T g are no longer the equilibrium ones, there is a priori a second contribution to AC p coming from the 
temperature dependence of the bottom of the well energy and of the vibrational properties of the TAP states 
for T > T g . In the harmonic approximation leading to Eq. these contributions come from the temperature 
evolution of en and of: 

K* = J dXp* (X)lnX, (21) 

where p*(X) is the density of eigenstates of the Hessian characterising the typical TAP states at temperature 
T. Since these TAP states are more and more stable as T decreases, dJC* /dT < 0. In the following, we will 
neglect these extra contributions and write AC p w T g dT, / dT\x as in the REM, but this particular point would 
need a more detailed discussion, since the change of K,* with temperature is in fact related to the temperature 
dependence of the high frequency shear modulus Goo(T), see section UlI El and Fig. [5l 

• The value of the exponent 6 and of the effective surface tens ion T p are not well known at this point. Based 
on exact calculations by Franz & Montanari in the Kac limit 105], one finds the naive result 9 — d — 1, which 
leads to I* cx Tn/TE. Recent numerical sim ulations ar e not incompatible with this result, albeit with a strong 
sub-leading correction to the surface energy jl06l . 1 1071 ] . Wolynes and collaborators, one the other hand, argue 
that 9 = d/2 ^ d — 1 and T = kT with k to a large extent independent of all molecular details. State to state 
fluctuations of T Q;7 must however play an important role, see below. 

• Because the surface to volume ratio goes to infinity when d — * oo, we expect To/f2<2 to diverge in that limit. 
Therefore I* also diverges in large dimensions, and the notion of "states" that span the whole system does make 
perfect sense in that limit. 

• The Adam-Gibbs argument also predicts the existence of a characteristic length Iag- The question they asked 
was: how large must a region be such that at least two metastable states can fit in? [3J] In other words, we 
want Y,{T)i d AG w In 2, leading to: 

In 2 ^ 1/d 



^ = lw)J (22) 

The above cavity melting argument only makes sense if Yj(T)i* d > In 2, corresponding to I* > £ag- Forgetting 
numerical constants (of order unity when d is finite), this inequality reads: 

E(T) < (l^y 6 , (23) 

which is surely satisfied close to Tr-, since E(T) — > 0. Note that even when the number of possible TAP states 
inside the cavity is large, the actual entropy of the cavity (as com puted from Z(R,T)) is small until R reaches 
£*, and only becomes E(T) for R > £*. Franz & Semerjian [To| quote E(R,T) = E(T)[1 - t/R] in the Kac 
limit (see below). 

• The cavity argument is consistent only if the number of TAP states inside the cavity is much less than the 
number of possible boundary conditions, otherwise these boundary conditions would not be able to select one 
given state. For R = £*, this inequality reads: 



e Sd 



> e n d e^(T) ^ I^n S ( T )«J-i-9 < Cdj (24) 



where Cd is a d-dependent numerical constant. This inequality is always satisfied close to Tk provided 9 < d— 1. 
So there is no contradiction in thinking that boundary conditions can fix the phase inside the cavity, even in 
short-ranged interacting models. 



D. Growth of amorphous order: order parameter and point-to-set correlations length 

The sections abov e makes clear that RFOT is a theory for an ideal glass transition where "amorphous long range 
order" emerges |223| . The glass phase should be characterised by an order parameter and a growing static correlation 
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length as the transition is approached. What do these concepts mean in the present context? The question of an 
order parameter for the glass transition is old. Back in 1983, Anderson (again) anticipated that there might be an 
hidden order parameter: Some - but not all - transitions to rigid, glass-like states, may entail a hidden, microscopic 
order parameter which is not a microscopic variable in any usual sense, and describes the rigidity of the system. This 
is the fundamental difficulty of the order-parameter concept: at no point can one be totally certain that one can really 
exclude a priori the appearance of some new hidden order. Insightful indeed: we now understand that within RFOT 
the order parameter is the overlap between two equilibrated configurations with the same boundary condition, and 
the correlation length is related to the "stiffness" of this overlap field. As outlined in Appendix A, the formal idea is 
to compute using replicas the thermodynamics of the cavity by constraining the overlap between the configuration of 
the system and the reference state a to be large outside the cavity. One can make analytical progress by considering 
a Kac model, where the range of interactions A becomes very large jl05l . ll09| . An interesting quantity is the resulting 
average overlap q(0; R) between the configuration at the centre of the cavity and the re ference state a , as a function of 
the radius of the cavity R. The quantity q(0;R) is called a "point-to-set" correlation [l04 Il05lllldl |. One expects on 
general ground that q(0; R) is large for small R and drops to zero at large R, since in that limit the system will explore 
configuration totally unrelated to a . The point-to-set correlation length £* characterises the extension of amorphous 
order. Montanari & Semerjian [llfjj | proved for a very broad class of systems that if the relaxation timescale diverges 
at a finite temperature then the point to set correlation length must diverge too. 

In fact, from our previous hand waving argument, we expect that the scale over which q(0; R) drops to zero is 
precisely £*. This is exactly what Franz & Montanari found in the Kac limit A — * oo; they indeed find a length 
I* beyond which the thermodynamical solution corresponding to a small q(Q; R) has a smaller free-energy than the 
solution corresponding to a high value of the overlap. This length f is found to behave as ATr-/(T — Tr-) when 
T -> T K , corresponding to Eq. (IBJ| with S(T -> T K ) cx T - T K , T (T K ) > and, in this Kac model, 6 = d- 1. 

Quite remarkably, another information naturally comes out of the calculation of Franz & Montanari, now about 
the finite dimensional physics close to Td [l05l Il08l |. The high overlap branch appears at temperatures above the 
dynamical MCT transition Td and remains metastable up to a length Id larger than I* , that diverges approaching Td 
from above as: 

'^ a G^J' "=i (25 » 

This means that typical TAP states, which should be unstable above Td, are stabilised by the presence of a frozen 
boundary condition, but only if the cavity is sufficiently small, i.e. R < id- For T — ► Td, the boundary condition is able 
to stabilise arbitrary large systems. One can understand this result in terms of delocalized soft modes: typical TAP 
states are indeed unstable above Td, but the corresponding soft modes have an extension ~ £d(T). This phenomenon is 
very similar to the one occurring in hard spheres packings close to the isostatic point (see [l 1 21 ] , where the importance 
of these soft modes is clearly discussed, and the analogy with the physics of glasses is underlined.) The existence 
of such a diverging length scale was in fact predicted within the contex t of MCT , on the basis of a very different 
dynamical calculation but with the same final prediction, Eq. (|2"5|) 42, 4~3l lll3llTT4T ]. Within MCT, Id is precisely the 
dynamical cor relat ion length £d, that measures how far a local perturbation affects the relaxation dynamics. As was 
argued in (42l . Ill5j . the value of v is expected to change whenever d < 8. An Harris- like argument on the fluctuations 
of the local free-energy suggests that if the transition survives, v should be > 2/d (see section IV CI for more about 
this). 

Finally, note that a distinct but possibly related static correlation length has been introduced by Kurchan and 
Levine [lll| . This is defined in terms of pattern repetition in a given configuration and, as the point-to-set length, 
should be infinite if there is amorphous long range order. 



E. Stability of the TAP states and high frequency shear modulus 



As the temperature is reduced, the typical TAP states not only have a smaller free-energy /* but are also more 
stable. The average curvature of the stability matrix is larger. Physically, this means that the fluctuations of the 
particles around their average position, (u 2 ) ~ T/G^T), is smaller not only because the tempera ture is smaller, 
but also because the stiffness, proportional to the high frequency shear modulus Goo(T), is larger. [224j Although 



TAP states become marginally unstable when T — > T d , in the sense that the smallest eigenvalue of the Hessian A n 



touches zero there, we expect that Goo in fact does not vanish at Td, 
the square-root singularity of the eigenvalue density: 



but has a square-root singularity inherited from 
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(26) 
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FIG. 5: High frequency modulus Goo (arbitrary units) as a function of e = 1 — T/Td for a semi-circle distribution of Hessian 
eigenvalues p*(A) with 2 shear modes per particles. Note that Goo increases substantially below T<j: this is a crucial ingredient 
in Dyre's shoving model (see section [VI A[) . We show on the same graph the vibrational contribution to the entropy /C* (up to 
an arbitrary vertical shift). 



This result can be established within MCT [117j | , and a replica calculation of a la Yoshino-Mezard (forbi dding 
inter-state transitions) should confirm this. Note however that there is no consensus on this point; for example, [1181 ] 
claim that G 00 (T ( ^) = in mean- field. 

This shear modulus gives a natural energy scale for the effective surface tension To, its singular behaviour around 
Td may turn out to be important for detailed comparison with experiments, see sections IIVBIVIA1 The fate of the 
square-root singularity in low dimensions and for realistic systems is however not known at this point . For elastic 
spheres close to the jamming point, for example, one finds that the shear modulus indeed vanishes jl!2| . 



F. Entropy driven melting: numerical evidence and fluctuations 

As we just discussed, the above cavity Gedanken experiment can in fact be used to make analytical progress. It 
can also be faithfully imple mente d in numerical simulations, for which the freezing of particles outside a spherical 
cavity is an easy operation. |225l ] One then monitors how the overlap q(0; R) between the reference state that acts 
as a frozen boundary condition and the configuration at the centre of the cavity varies with R. The point-to-set 
correlation length £* can be defined, for example, such as q(0; R > £*)/q(0; a) < 0.5. This length scale is clearly found 
to increase substantially when temperat ure is red uced below X^, at least in the temperature range where numerical 
simulations can be equilibrated [l!9lll20| (see also jl2l| for a related effect). Quite interestingly, the shape of q(0; R) as 
a function of R also evolves very significantly from a simple exponential decay at hi gh t emperatures to a "compressed 
exponential", ln<?(0;i?) oc — with £ > 1 at lower temperatures. As argued in jl20l |. the compressed exponential 
shape may be a signature of fluctuations. One type of fluctuation that we already mentioned comes from the effective 
surface tension between different TAP states, T a;7 . For a fixed external state a, it is indeed conceivable that for a 
given R some 7 states are good matches and have a particularly low T Q;7 . 



1. Surface energy fluctuations and renormalisation of 9 

To understand how these fluctuations might strongly affect the physics, we now write the partition function for the 
mobile cavity surrounded by the pinning state a as: [120] 



Z a (R,T) = e 



- e -/3R d f a 



7^a 



(27) 
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The probability pa to leave state a is therefore: 



-f3B d ^-f3B a T a; ,, 

Introducing Af a {f, T) = J(/ — / 7 )<5(T — T Q;7 ), the sum can be written 

£ e -/3fl%-/« e T QT = J df J dTe -W d f-W^M a U,T), (29) 

= JdfJdre-p Rd f-P R<>T+Rd ^^p a (r\f), (30) 

where in the last equality we have defined p Q (T|/) = Af a (f, T)/7V(/), and A/"(/) = exp[<r(/, T)] is the number of states 
with free energy /. P a (T\f) is the fraction of states 7 (inside the cavity) with free energy / and effective interface 
tension T. The integrals above can be simplified using the saddle point method, which is a very good approximation 
even for rather small values of R. The integral over / picks up the most probable value / = /*, as above. Because 
of the negative exponential term in R , the integral over T is dominated by the lowest values of T allowed by the 
distribution P a (T\f*). There are two possible cases: (a) If this function vanishes for T < To (a), one finds that 
up to sub- leading terms one can forget about the fluctuations of T Q:7 and use To(a)R e in the cavity argument. 
However, if (b) there are arbitrarily small effective tensions T, the value of 9 can be renormalised. For instance in 
the case where P a (T\f*) ~ exp[— (T (a)/T) y ], where y is a positive exponent, one that the effective surface energy 
is now proportional to R e , with 9' = y0/(y + 1) < 6, and again an a dependent effective value of To. This is quite 
important because the mosaic length I* is determined by this renormalised value of 9. Wolynes et al. have obtained 
9 = d/2 < d — Ion the basis of a "wetting" mechanism that softens the interface. It is unclear to us whether this 
wetting mechanism and the existence of rare, but arbitrary small surface tensions, in fact describe the same physical 
effect (in which case, one should argue why y = d/(d — 2)). 

We have shown how the fluctuations of T a;7 for a fixed a can be accounted for. But it may still be that the effective 
surface tension To (a) strongly depends on a - for example, a could be exceptionally stable locally, corresponding to 
a large To (a). The average overlap q(0; R) receives contributions from samples that are locked in state a, i.e.: 

? (0; R) « / ^oP(To) 1 + exp(s(r) ^_ Tojyvr) (3D 

where P(T ) is the pro babil ity (over a) to obtain a renormalised surface tension T , which may depend both on 
R and T. As shown in |l20l |. this assumption can be used rationalise the compressed expone ntial shape of q(0:R) 
found numerically. In a recent work P( Tp) h as been directly measured using forced excitations [l07j . The results are 
compatible with the shape surmised in [12C| and underline the importance of the fluctuations of Tq. 



2. Free-energy fluctuations and "locally preferred structures" 

A second, very important source of fluctuations comes from the local free-energy of the selected TAP state. In 
the limit where the size of the cavity £* goes to infinity (i.e. close to Tk or in the mean- field or Kac limit), only 
states with free-energy density /* have a significant probability of being observed. But when I* is finite, there is some 
probability to observe exceptionally low free-energy states; when f — f* is small, the distribution is Gaussian with a 
width cx l*- d / 2 . Let us assume for simplicity that the complexity <r(f, T) takes the following form (by analogy with 
Eq.© above): 

a(f, T) «/?*(/- /min(T)) - y (/ - Jmin(T)) 2 + . . . f> / min , B > 0, (32) 

with / m in(T) = e m j n + f v a(T), where f V ib(T) is a certain function that describes the free-energy of the (possibly 
an-harmonic) vibrations. This explicit choice leads to: 

.r=/ mi n(T) + ^^ + ..., E(T) = fc^ + ... (33) 

and / — /* = u/\/ B£* d , where u is a Gaussian variable with unit variance. As noticed above, the jump of specific 
heat at T g is related to B as: AC P = ftg/B, if we neglect the vibrational contribution (see discussion in section ITlI C[) . 
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A TAP state a with a lower free energy (u < 0) is more stable against fragmentation, not only because the typical 
states that can 'nucleate' now have a lower free-energy / < /*, but also because the effective surface tension is 
expected to increase with the stability of the state, and therefore with /* — / (see section IIIIE|) . This means that the 
size of the cavity beyond which this particular TAP state becomes unstable is increased. The balance between the 
free-energy gain and the surface energy loss now gives, to first order in u: 



t 



Au 2(3f3 g 



(34) 



where A > is a numerical constant. Using £* = (/3To/E(T)) 1 / d_e , the correction term is found to be small provided 

I3 2 K -P 2 \ 29 ~ d „ (/3T ) d 



2flSg J <<: ^Cf' (35) 

which is valid arbitrary close to Tjc only if 9 > d/2, and marginally valid when 9 = d/2 which happens to be KTW's 
conjecture. If 9 < d/2, free-energy fluctuations in fact dominate the surface energy cost, and quite interestingly, leads 
to an effective value of 9 stuck at d/2. That 9 — d/2 plays a special role was in fact argued early on by KTW; it 
corresponds to £* diverging as (T — Tr-) -2 ^ close to Tk, as expected for a disorder dominated phase transition Q. If 
indeed 9 = d/2, the criterion Eq. I|35[) is never satisfied unless AC P <C 1, which is not the case for fragile glass formers 
(see Fig. 3-b and below). We expect in this case fluctuation effects to be strong, as one may have anticipated from 
the "wetting" argument used by KTW to obtain 9 = d/2, which also predicts that the width of the domain walls is 
itself £* , i.e. that domain walls are in fact everywhere! 

When £* is finite, there is an exponentially small probability to even observe the lowest possible free-energy states, 
with / = / m in(T). These states can be called "locally preferred structures" (LPS) by analogy with theory of glasses 
based on the existence of locally preferred pack ings (for a review, see [l2j and section [VI Bl and for recent interesting 
ideas about the nature of these LPS, see |l22j ). These maximally stable states are robust against fragmentation up 
to size: 

t max (T)=t(T)i^^j , (36) 

which becomes much larger than £* close to Tk ■ The probability to observe such a large droplet is however exponen- 
tially small, and given by: 

P(Cax) « exp[-CL£(r)] (37) 
These fluctuations effects should have important consequences for the dynamics, in particular for fragile systems 
(because the quantity in the exponential behaves, for a fixed T/Tk, as ACp S ^ d ~ 8 ^). We will come back to this point 
in section ITV CI 

G. Summary 

The point of the above section was to establish, on the basis of purely thermodynamic arguments, that the notion 
of stable TAP states with a finite configurational entropy £ makes sense physically only if one zooms on a region 
of space not larger than a certain length £*, that diverges in mean-field or when £ — > 0. Above £* , the system is a 
liquid, i.e. thermodynamically micro-phase separated on the scale £* into TAP states. Below £*, on the other hand, 
the system is an ideal glass even when T > Tk- 

In the course of the argument, one has to introduce the somewhat hazy notion of a surface energy between two 
different TAP states, and assume that it grows as TqR . There is little agreement at this stage about the value of To 
and 9. Mea n-fie ld calculations and numerical simulations seem to favour 9 = d — 1, possibly with strong finite size 



corrections [107j . whereas Wolynes et al. recommend 9 = d/2. We have seen that fluctuation effects, absent in mean 



field, could renormalize the value of 9, and that 9 = d/2 plays a special role in that respect. 



IV. DYNAMICS IN THE MOSAIC STATE 
A. Decorrelation and activation 



Let us re-interpret the above results from a dynamical point of view. We have seen that if the configuration outside 
the cavity is frozen in state a, the particles inside the cavity have an overwhelming probability to stay in state a 
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whenever R < £*. Dynamically, some low energy excitations may sometimes be explored, but the system always 
reverts back to the a configuration. The dynamical structure factor C{i) does not decay to zero. The relaxation 
modes of length R < £* cannot be used to restore ergodicity in the system, and it is therefore self-consistent to 
assume that for these length scales, the environment is frozen, or more precisely fluctuates around a well defined 
configuration. The motion on these length scales corresponds to a generalised 'cage' effect and contributes to the 
so-called (3- relaxation, i.e. the approach to the plateau value C(t) as q* . Interestingly, the back and forth motion 
of groups of mo lecu les between two or more configurations has been clearly observed using mesoscopic investigation 
techniques, see [l23j . 

When R > £*, on the other hand, the configurational entropy becomes so large that the system ends up in a 
completely different state 7 with only an exponentially small probability to come back to state a, even if the outside 
of the cavity is, by fiat, frozen. Therefore, the relaxation time of the system is the time needed for a cavity of size ~ £* 
to leave state a - once this is achieved, the configuration inside the cavity will most probably never come back. The 
dynamical structure factor now decays to zero, and the assumption of a frozen environment is no longer self-consistent 

A precise characterisation of the dynamical processes that allow the particles in the cavity to leave state a is an 
open problem. These processes must however satisfy some constraints. First, they have to be compatible with time 
reversal symmetry. Second, the transition matrix should be such that many arrival states 7 can be reached once the 
system has left a, otherwise the system would not gain configurational entropy. Third, the dynamical process has to 
be compatible with the fact that TAP states are locally stable. As discussed above, a crucial point is that on scales 
R less than I* the system is like an ideal glass, and therefore the barrier to escape from the state a should grow with 
R. On the ba sis of these remar ks, a reasonab l e gu ess (inspired from known results in the context of pinned systems 
[l24i 1 125L Il26l | and spin-glasses jl27l Il32l Il33l [l34j |) is that the arrival state will be reached through barrier crossing 
and the lowest energy barrier to significantly alter the configuration in a cavity of size R grows with R as a power-law: 

A(R) = A i?^, (38) 

where Ao is an energy scale that presumably depends on the departure state (again, lower free-energy states are more 
stable and should have a larger Ao) and ip another exponent. A positive entropic correction is possible. This would 
be due to the fact that in order to leave the state a the system not only has to jump a barrier but also to find the 
target states between all possible configurations. The higher the temperature the larger the number of the target 
states relative to all available states. As a consequence this entropic contribution to the barrier should decrease at 
fixed R when increasing T. 

Naively, one may think that the reconfiguration events proceed by nucleating a droplet of another phase, leading 
to tp = 6 and Ao = To- This was KTW's original surmise, and is still advocated in Wolynes' recent papers (see [ll| 
for a review). More generally, one expects that the barrier should be at least equal to the excess energy of the final 



state, leading to ip > 6, or ip = 9 and A > T (127j |. Needless to say, a prec ise analyt ical or numerical determination 
of Ao and ip is currently beyond anyone's ability. Recent numerical results |l06l . Il07| suggest ip as 1 and 9 as 2. This 
is actually quite puzzling in view of the expected boun d 9 < . ip. At this stage it is unclear whether the violation of 
the bound 9 < ip is related to the way 9 is measured in jl06l 1 1 07t j or signals that something is missing (or wrong) in 
the theoretical description. 

Even the nature of the activated events that determines the order of magnitude of Ao is quite obscure. Is Ao 
proportional to T, as assumed by Wolynes and associates, or is Ao more akin to an elastic energy, and therefore 
proportional to Goo(T)? These uncertainties are related to our lack of understanding of the physical nature of the 
activated events. Should one picture them as the activated fluctuations of the grain-boundaries between glassites that 
sweep larger and larger regions of the bulk until £* is reached? Or is it the nucleation of small droplets within the bulk 
that grow and percolate? Or are these two pictures in fact equivalent because, as implied by the "wetting" argument 
of KTW that leads to — d/2, the width of the interfaces are themselves o f ord er 1*1 What is the role played by 
"excitation chains", of the type proposed by Glotzer et al. [l28j and Langer [l29j |? 

In any case, the typical time needed to decorrelate a cavity of size £* is obtained from Eq. (|3"8|) as: 

/ A £*^ \ 

T a ~ r(t) = to exp ( —5- — ] , (39) 



T 

where we neglect all fluctuations at this stage (see section llV Cj) . as well as a possible power-law prefactor £* z in 
the Arrhenius formula. From the above arguments, it is clear that the a-relaxation time of the system is t(£*): 
smaller length scales are faster (energy barriers are lower) but unable to decorrelate, whereas larger scales are orders 
of magnitude slower so that the evolution on these scales will be short-circuited by a relaxation in parallel of smaller 
blob s of size £* . The divergence of r a when £* — > 00 is completely in line with the analysis of Montanari & Semerjian 

DjU- 
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It is interesting to reformulate the above argument as a discussion of finite size effects. A cavity of size R < £* should 
relax to equilibrium in a time t(R) -C t(£*), although this relaxation is incomplete in the sense that C(t ^> t(R)) — > 
q*(R) > 0. As the size of the cavity increases beyond £*, the relaxation time should increase and then saturate around 
t(£*), whilst q*(R) — > 0. This scenario is expected for temperatures T < Td, and preliminary numerical results 
seem to confirm it |l 35l j . Certainly, this is a major point that needs to be carefully tested in the future numerical 
investigations. When T > Td, on the other hand, more complicated size effects should be present. One expects that 
the cavity equilibrates through non-activated relaxa tion channels corresponding to collective unstable modes of size 
£d ~ (T — Td)~ u (see discussion in sectio n IIIID I and [105[). If R < £d these channels are not available and relaxation 
must be slower, not faster (see [l30L Il3l| ). Thus one expects that the convergence towards the bulk relaxation time 
r a is from above for T > Td and from below when T < Td- However, the situation might be more complicated still 
because of the presence of activated events above T d (see below and section IVCj) . 

Before moving on, it is worth addressing a worry expressed by Langer: is it consistent to take into account the 
exponentially large number of TAP states to compute a static £* and then to claim that after a time t(£*) needed 
to visit only a few of them the whole structure has evolved? Is the full configurational entropy S(T) meaningful 
physically, or only a small fraction of it? This is however a standard problem in statistical mechanics: although no 
physical system ever visits the exponential number exp(iVS) of accessible states, we know that after a relaxation time 
much smaller than tq exp(A^S), thermodynamical equilibrium is reached and the use of the full entropy is warranted. 
What is important is not the number of states that are effectively visited, but rather the number of states that can 
be visited in the course of dynamics. This discussion would however be relevant to determine the correct power-law 
prefactor £* z in the above Arrhenius formula for the relaxation time. 



B. Consequences: Adam-Gibbs and Vogel-Fulcher 

Inserting the expression of £* given in Eq. dTHJ) into Eq. we finally obtain: 



T T \ TS(T) 

where Cd is a numerical constant which is difficult to estimate based on the above hand-waving arguments. The 
important point, though, is that the Adam-Gibbs inverse relation between relaxation time and configurational entropy 
appears very naturally. Physically, it reflects that configurational entropy is the factor that limits the growth of 'hard' 
glassy regions (glassites) which resist shear. Since these regions must evolve for the system to flow, and since the 
energy barrier must grow with the size of the glassites, the Adam-Gibbs correlation follows. 

The KTW conjecture ip — 8 = d/2 allows one to recover precisely the Adam-Gibbs relation, whereas the naive 
values ip = 9 = d — 1 lead to a stronger dependence, as lnr Q ~ E~ 2 in d = 3. The values of the prefactors A , T 
(and their possible fluctuations) are obviously relevant for a more quantitative comparison with experimental data. 
Acco rding to Wolynes et al., Ao = To = kT, with k a numerical constant nearly independent of all molecular details. 
[226| Using the REM like model above for the configurational entropy then leads to a modified Vogel-Fulcher form 
for the relaxation time: 

Tq 2C d n 2 p 2 g C d n 2 T K 

ln To= ACM-P 2 ) ~ AC P (T-T K ) ( 41 ) 

The calculation of the fragility parameter to, defined in Eq. Jl} is quite interesting. Within the Wolynes choice of 
parameters one finds, independently of the precise shape of S(T): 



m = _ T 91og 10 r Q 



dT 



mo w (42) 



where mo = 16 is the number of decades separating the microscopic time scale and the conventional relaxation 
time defining a glass. But within the same Wolynes framework S(T ff ) is nearly universal, since by definition of T g , 
too = C'dK 2 /Ti(T g ) (with too = In 10 too) The final result is therefore a simple relation between the fragility m and the 
jump of specific heat ACL, first derived by Xia and Wolynes Q: 



pi 

m m Q 



CdK 



AC P . (43) 



Experimentally, the slope between to a nd A C P is of order 7 (see Fig. 3-b), leading to S(T g ) sa 2ks per particle, 
independently of the values of n and Cd- [227] This result is quite reasonable indeed: as mentioned in the introduction, 
molecular glasses are characterised by an excess entropy of a few Ub at T g (see Fig. 1). 
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Another possible choice would be Ao = To = nT g , independently of temperature. The fragility parameter is now 
found to be: 

m = 2m Q + ^^AC P , (44) 

i.e. a non zero intercept in the relation between m and AC p can be accommodated by relaxing the strict proportionality 
between Ao, To and temperature. 

Within both hypothesis, the relevant mosaic length scale at T g is also universal, and given by: 

e ^-iMT (45) 

corresponding to a number of particles inside the 'critical' cavity N*(T g ) — (4tt/3)£* 3 w 7OC3 in d = 3. Taking C3 ~ 1 
for lack of other natural choices, one finds N* ~ 70; less than a hundred particles have to move together to make the 
system flow at the gla ss temperature. Th is is not inconsistent with various experimental estimates of this number in 
molecular glasses [4t| Il38l Il39l Il40l Il4l| . Note however that the associated length scale at T g is still very modest, 

Both the above choice for Ao and To is however subject to quibble and is in our opinion an open and pressing 
problem. That the interface energy between two TAP states, and the energy barrier between them, decrea se wh en T 
is reduced, is not intuitive at all, in particular in view of the temperature dependence of the shear modulus. 228] Why 
should the microscopic mismatch between two TAP states decrease at low temperature wher eas the sh ear modulus (a 
proxy for a local barrier when multiplied by a microsco pic volume a 3 ) does just the opposite? [I44lll45| This discussion 
is in fact related to Dyre's "shoving model" jl46l . Il47| |. where the super- Arrhenius behaviour of the a-relaxation time 
is entirely attributed to the growth of the high frequency shear modulus, through: lnr Q /ro oc Goc(T)/T, with no 
contribution of the configurational entropy. We will come back to this point in section IVI Al 

C. Dynamic fluctuations: stretched exponentials and facilitation 

As we explained in section UlI Fl we expect strong local fluctuations to be present in glassy, amorphous systems. 
These fluctuations affect both the size of the glassites £* and the energy barrier parameter Ao . The relative fluctuations 
of the local relaxation time are given by: 



= In — x 



SA 81* 
A^ + ^V 



(46) 



Particularly low free energy states have a large £* and presumably a large A ; however T and A also fluctuate on 
their own, with relative fluctuations of order £*~ UJ y where w is another unknown exponent. The free-energy induced 
fluctuations, on the other hand, are of order T 9V /AQ;r- d / 2 . The total variance of the local fluctuations of relaxation 
time is therefore given by: 




Av£*~ 2u + At Tg ^ Cp t~ d 
+ yl/ T 2 £(T) 2 



(47) 



where Af and Af are numerical constants with subscripts tagging their physical origin. As is well known, local 
fluctuations in the relaxation time induce both a stretching of th e relaxation function, and a decoupling between 
viscosity and self diffusion, often called SER violations [23L l40l . fl48l ] . A rough, but useful interpolation formula r elate s 
the exponent (3 of the stretched exponential relaxation to the width of the relaxation time distribution through [l49j ] : 

'1/2 

(48) 

This relation is interesting. If we believe that £*{T g ) and S(T S ) are indeed universal, then very fragile liquids with 
AC P > E(T g ) should be such that, in d = 3: 




« T »»~ s <W;diW~V^ <49) 
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where we have used the linear relation between m and AC P obtained above. This dependence appears to be 

compatible with the compilation of results provided in where the proportionalit y con stant is found to be « 3, 
leading to a Af ss C3/9. An inverse dependence between m and (3 is also reported in [l5Cj . 

In the same limit of fragile liquids where the frcc-cncrgy fluctuations should dominate, one finds that the stretching 
exponent (3 decreases when temperature is decreased, as seen experimentally. More precisely the relation reads 
/3(T) oc Tl* 6 / 2 -^ in the limit /3{T) < 1. In particular, /3(T -> T+) -> whenever ip > 6/2, in particular if tp = 9. 

The typical fluctuations of the relaxation time are adequate to understand the shape of the correlation function 
around its inflexion point. But the long time behaviour of the relaxation will be domi nated by very rare, but extremely 
stable states. One can come up with a droplet argument, inspired from spin-glasses |l5l| . for the relaxation function 
that illustrates this point. Let us consider droplets of all sizes £ around a given point in space. Each droplet can have 
a free-energy / not necessarily equal to /*. For a given £ and /, the probability to be still in the initial configuration 
after time t is: 

_ t 

PiW)= 1 + *i</^\ */(/) « exp [(E(T) - /?(/* - f))i d - PT f] , (50) 

where t{1) = r exp(/3A ^), and we neglect the fluctuations of To and Ao for simplicity, which is a better approxi- 
mation for fragile materials according to the discussion above. The above equation tells us that after a time ~ t(£), 
the probability to be still in the initial state reaches its equilibrium value (1 + z^) -1 , where zi is the weight of all 
other droplets of size £ that can nucleate. In principle, this quantity should read: 

Mf) = I df'ffW'-M'-M'-P***, (51) 



and a saddle point approximation leads to the above expression with /' w /* and cr(/', T) w S(T). 

The correlation function is the probability that no droplet, whatever its size, was able to change state. Therefore: 



at) « n 



(a(/,T)-/3/)f d 

df Z(£) P ^ 



(52) 



This expression is complicated to analyse in general but is quite interesting since it should in principle provide an 
approximate description of C(t) in the (3 regime C(t) sw q* in terms of droplet excitations. For very large lnf, a saddle 
point estimation leads to the following result: 

/Tin A 

lnC(t) ~mt»i -AC p T 2 (f3 K - f3f J , (53) 

i.e. a relaxation slower than any stretched exponential, that is dominated by the extremely low free-energy, "preferred 
structure" states with / = / m i n - Note that: 

• This asymptotic decay law allows one to define a characteristic relaxation time t^. From C(t > Too) < £^ one 
finds: 

lnr^cx^-TK)- 2 ^, (54) 
which diverges with the same exponent as lnr Q whenever 9 = d/2, independently of the value of ip. 

• This slow relaxation tail is expected to be present in the range Td < T < Tq. Above To, even the deepest states 
become unstable (see also section IV CI below) . 

Note that when an activated event takes place within a droplet of size I* , the boundary conditions of the nearby 
droplet changes. There is a substantial probability that this triggers, or facilitates, and activated event there as 
well, possibly inducing an "avalanche" process that extends over the dynamic correlation length scale £<j > £*. The 
dynamics on length scales < £* is, within RFOT, inherently cooperative, but one may attempt to construct a coarse- 
grained description of the dynamics beyond £*. What is the appropriate model? It seems to us that even if some 
facilitation mechanism is indeed highly plausible, the coarse-grained model should be built as kind of contagion model 
of activity with Pois son l i ke ac tivat ed in i tiato rs, rather that a Kinetically Constrained Model with strictly conserved 
mobility defects [36l Il55l . Il56l | (see [l57l Il58| for a recent discussion of this point, and section |VI CI below') . In any 
case, the relation between the dynamic correlation length defined for example through three- or fo ur-po int point 
correlation functions 0, and the mosaic length £* is at this stage an important open problem (see [l39l Il4p| ). 
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D. Two new predictions of RFOT 

1. Energy relaxation 

The mosaic picture of RFOT theory suggests an interesting scenario for the relaxation of the energy per particle 
after a sudden quench from high temperatures, T 3> ?d, which in turn may provide a test for the theory and a way to 
measure (numerically and experimentally) important quantities such as Tq(T) and the exponents 9 and tp. Although 
some predictions of RFOT in such an out of equilibrium aging regime were made in fTTI j , we believe that the discussion 
below is new. 

The idea is that after a sudden quench, TAP states appear locally but initially have a very small extension £. After 
a rather short stage after the quench where the initial excess energy disappears, a slower regime where the TAP 
droplets have to grow to reach their equilibrium size £*. The residual excess energy is primarily concentrated in the 
interface regions. If the typical size of the droplets at time t w after the quench is £(t w ) <C £*, it is reasonable to expect 
that the excess energy e per unit volume is given by: [l32j | 

e (^,T)- e *(r)a To ™j; )]g (55) 

where e*(T) is the equilibrium energy at temperature T. If we assume that the droplets have to cross the same kind 
of energy barriers to grow as the ones that are dynamically relevant in equilibrium, the time dependent droplet size 
should be given by: 

When £(t w ) gets close to £*, the final relaxation stages correspond to smaller scale rearrangements and are expected 
to be much faster. Therefore, a natural conjecture for the energy relaxation of the mosaic state is: 

with T{x — > 0) equal to a numerical constant, whereas T(x — » oo) decays quickly to zero. The decay of the excess 
energy should therefore be an inverse power of the log of time, at variance with models of conserved mobility defects 
where the decay should be faster: ~ t w d ^ z . The claim is that the scaling form Eq. ([57| should account for different 
temperatures and allow one, if correct, to extract useful information on Tn( T) an d Aq(T). 

Preliminary numerical simulations seem to be compatible with Eq. (|57[) [159]. Careful calorimctric experiments 
might allow one to test this prediction for real supercooled liquids. However, a major technical obstacle is to be able 
to reach very fast cooli ng ra tes, or else one should generalise Eq. (j5"Tj) to an arbitrary temperature cooling scheme 
T(t w ), as was done in [160j for a similar problem. For example, one could study regular quenches with constant 
cooling rate F. Up to logarithmic accuracy, the system is equilibrated up to a length £(T) given by: 

/T-lnJ-V'* 

^(a^t) • < 58 > 

where T* is the crossover temperature at which the energy barriers start growing (in fact, a slightly more precise 
formula would replace T* / Aq(T*) by Ti nt / A$ (Ti nt ) where Ti nt (T) is some intermediate T-dependent temperature 
between T* and the final temperature T.) As a consequence the residual energy should scale as: 

6-d 

e(T,T) - e*(T) oc T„(T) ( j (59) 

If one believes that most of the excess contribution to the age dependent volume or dielectric constant also comes 
from the interfaces, then similar scaling laws should hold for these quantities as well, with 8 replaced by d s , the 
fractal dimension of the interfaces. If d s = d— 1, the decay of these quantities as (lni^,) -1 /^ sh ould give direct access 



to the value of ip. These ideas could be used to interpret the volum e experiments of Kovacs [16l| or the dielectric 
experiments of Leheny & Nagcl [162] or of Luckenheimer et al. |l63l ] 
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FIG. 6: Sketch of the shear dependent inverse viscosity n' 1 = 7/cr as a function of a. For a < a*, 7/cr w 1/77;, the zero-shear 
(linear) velocity. For a > a*, some significant shear thinning is expected. If the RFOT scenario is correct, the curves for 
different temperatures Ti < T2 should look as drawn (at least for low enough temperatures when Goo(T) is approximately 
constant), with lower temperatures corresponding to larger viscosities and lower shear thinning stress. Other theories lead to 
the opposite behaviour, where a* increases upon lowering the temperature. 



2. Non linear flow curves 



Another inte restin g consequence of the putative mosaic structure of the liquid is the non-linear response to an 
external shear. |229l | The elastic energy stored in a volume R d due to a shear a is given by cr 2 i? d /2G 00 . If the elastic 
energy of each glassite is smaller than the typical energy barriers Ao£*^ that the liquid has to spontaneously cross to 
flow, the rheology will be Newtonian, with a shear rate 7 = cr/rji, where r/i = Goo^a is the linear viscosity. But if this 
elastic energy becomes larger than Ao^*^, the external shear lowers the barriers to flow and speeds up the flow. In 
this regime, the highest barriers correspond to a smaller droplet length scale i a such that both energies balance: 

1 

Ao# ~ a 2 ^a 2 /2G oc - = (^^) '~* , (60) 
such that 77 becomes strongly reduced by a (shear-thinning), as: 

r)(a < a*) w m ; r)(a > a*) oc rjf M , with Y(a) = (^-j < 1, (61) 

where a* is the crossover shear stress where the rheology becomes non- linear, corresponding to i a — £*. Assuming 
for simplicity ip = 0, and To ~ Ao ~ Goo, one finds: 

ct*(t)~ v/rs(r)Goo. (62) 

This predicts, rather non-trivially, that the cross-over shear stress where the viscosity starts plummeting decreases as 
the temperature decreases and the viscosity itself increases, at least if Goo is approximately constant. As T — > Tr-, 
<t*(T) ~ yJT — Tk ■ The shape of the flow curve should schematically look as in Fig. [6l at least at low enough 
temperatures where the variation of Goo(T) with temperature can be neglected. But since in many fragile systems 
Goo(T) increases quite substantially between T* and T g . The behaviour of a* with temperature and the corresponding 
flow curves could then be quite different. 

Intuitively, the reduction of a* as temperature is decreased is intimately linked to the growth of the amorphous order 
correlation length £*, below which the system behaves as a true thermodynamic glass. Experimental confirmation of 
such an evolution with temperature would be a very valuable confirmation of the premises of RFOT, in particular 
because alternative pictures lead to the opposite evolution with temperature. 
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E. Summary 

RFOT asserts that the fast growth of the a relaxation time is primarily due to activated rearrangements of glassy 
droplets of ever increasing size £* , in turn induced by the reduction of configurational entropy. This framework allows 
to account naturally for the Adam-Gibbs correlation between dynamics and excess entropy. It allows one to rationalise 
the relation between fragility, relaxation time broadening and specific heat jump. The existence of a non-trivial glassy 
length scale £* should show up in the non-linear rheological properties of supercooled liquids, that should become 
more pronounced as £* increases, i.e. at lower temperatures. We also expect that the growth of the mosaic length 
£(t w ) in out-of-equilibrium situations should have interesting observable consequences. 

V. THE MYSTERIOUS MCT/RFOT CROSSOVER 

As we repeatedly argued above, one of the very strong selling point of RFOT is that Mode-Coupling Theory is 
naturally embedded within the theory and describe the initial stages of slowdown of the dynamics, at temperatures 
higher than the dynamical (MCT) transition Td- MCT offers several precise predictions for the scaling of the relaxation 
function, and for the divergence of the relaxation time (see [26], [U HH and section III B 31 for a short summary). 
Within an energy landscape interpretation, at high temperature the system is typically near unstable saddle points. 
The slowdown is due to the fact that these saddles become less and less unstable as T approaches Tj, from above, 
see (63l . Il65| |. As shown in (42j, relaxation through these unstable modes involves a growing number of particles: the 
dynamic correlation length grows like ~ (T — Tc)^ 1 ^ in mean-field. 

In finite dimensions, however, the relaxation time does not diverge because at some point activated events become 
more efficient a relaxation channel than unstable modes. This is precisely in line with Goldstein's early insight 
about the existence of a crossover temperature T* separating "free flow" for T > T* from activated dynamics for 
T < T*. But the details of the crossover are still poorly understood theoretically, and potentially very constraining for 
RFOT. In spit e of several attempts, only very qualitative RFOT predictions for this dynamical crossover are available 
1 1671 Il68j | , see also [28j . In the following we will review the evidence, or lack thereof, of a clear dynamical 
crossover around T* of the kind envisaged by RFOT-MCT. We will also propose some very crude ideas to rationalise 
a somewhat puzzling situation. 

A. Encouraging items 

Many hints of a crosso ver i n the dynamics around a temperature such that r Q ~ 10~ 7 sec. have indeed been 
reported in the literature [l69j . We show for example in Fig. [7J the relaxati on t ime r Q as a function of temperature 
in a representation where the Vogel-Fulcher law becomes a straight line [l70l |. This clearly reveals a change of 
behaviour around r Q ~ 10 -7 sec, which is even more apparent in the bottom panel of Fig. [JJ where the ratio of r a 
to its low temperature Vogel-Fulcher fit is represented. Several authors have insisted that the effective energy barrier 
A(T) = TlnTa/To starts growing significantly below a certain temperature, that we identify with T* , as: 

A(T) = A > + Ao(l-^y, (T < T*), (63) 

where A> is the temperature independent barrier above T*, and tp « 8/3 0, HI or ip = 2 [H - see section ECE 

The literature lending experimental or nu meri cal support to MCT is abundant and reviewing it is far beyond the 
scope of this paper; we refer to [26|, [2?], [H, Il7l| for overviews. We note that MCT seems to account well for some 
aspects of the initial slowing down of the dynamics, in particular the appearance of a two step relaxation decay with 
a non trivial p regime that abides Time- Temperature superposition (see Eq. (|13p . The phenomenon of dynamical 
heterogeneity, that experim entally or numerically emerges already within the MCT regime, is also correctly captured 
by the theory [42|, l43l Ill4j - at least qualitatively. This is actually a very stringent test for MCT; as a matter of 
fact the founding fathers of MCT never believed that the theory could predict the growth of such a dynamic length 
scale! But had that not been the case, it would have been hard to argue that MCT is a physically sound theory for 
supercooled liquids. 

However, quite unfortunately, none of the quantitative MCT predictions can be tested beyond any quibble. This is 
because all these predictions only become exact (a) in mean field and (b) extremely close to the dynamic singularity T c i 
(see ref. (96J). Not only critical fluct uatio ns are expected to renormalize all MCT predictions whenever the dimension 
of space is less than d c = 8 [12, 0, Ill5| , but the dynamic singularity itself in fact disappears and becomes a mere 
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FIG. 7: Evidence of a crossover temperature T* such that r Q ~ 10 -7 , in a representation where Vogel-Fulcher law becomes a 
straight line. Lowe r pan el: same data, but divided by the Vogel-Fulcher fit adequate for T < T* . Here the crossover is most 
conspicuous. From |170( ]. 



crossover. Therefore MCT cannot claim victory before we have a full understanding of these two issues. The details 
of the MCT-RFOT crossover is clearly part of this predicament. 

There has also been some numerical effor ts to investigate directly the nature of the saddle points of the energy 
landscape that are probed by the dynamics [165]. The results are compatible with the idea of an energy threshold 
(or more probably an energy crossover) above which saddles are indeed mostly unstable. The average index (i.e. the 
fraction of negative eigenvalues of the Hessian matrix) of the saddles as a function of the energy per particle is shown 
in Fig. [5J both for realistic models of glass formers and for the p-spin model. As we mentioned above, the vanishing 
of the saddles' index is associated to the divergence of a correlation length, which in turn should be associated to a 
delocalisation of th e cor responding eigenvector of the Hessian matrix. Although direct evidence for this is lacking (to 
say the least - see |l72j |). numerical investigations of dynami c cor r elations within the expected MCT region indeed 
reveal the existence of a growing dynamic length & — £d{T) jl3ll . Il39l . Il73l . Il74 Il75l |. The quantitative agreement 
with the predictions of MCT is however still under debate. From that perspective, one important qualitative result 
is that dynamical correlation length appears to increase quite fast (as a power l aw of rel axation time) in the MCT 
region, crossing over to a much slower, logarithmic growth at low temperatures jl39l Il74| . as expected if activated 
dynamics sets in. Thus one indeed finds some indication of a cross-over in the structure of dynamical heterogeneities 
compatible with the MCT-RFOT scenario. 

In summary, MCT certainly reproduces qualitatively many of the physical phenomena pertaining to the dynamics 
of moderately supercooled molecular glass formers. However, it cannot claim victory on the basis of its quantitative 
predictions: some of them are indeed remarkable but there are also notable failures that will be discussed in the 
following section. One could argue that this is inherent to the fact that in finite dimensions the MCT transition must 
necessarily morph into a cross-over. This is indeed reasonable, but it also implies that validating the theory by testing 
its quantitative predictions is a hopeless task. We believe that one should instead devise demanding tests to ascertain 
that the physical phenomena taking place around Goldstein's cross-over temperature T* are indeed the ones encoded, 
even in a oversimplified way, by the MCT formalism. A truly smoking gun evidence should result from the study of 
finite size effects. These should be strongly anomalous, in the sense that above T* , smaller system should relax more 
slowly than larger systems, in particular when R < t&, since unstable modes should be stabilised by the boundary 
conditions (see section fill D|) . Instead, below T*, the relaxation mechanism should consist in activated dynamics and 
within RFOT one expects that smaller systems should rel ax faster, at le ast until R « £* . Some evidence for such an 
anomalous size dependence has indeed been reported in [l30l Il3ll Il35j | , but more work in that direction is clearly 
needed. 
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FIG. 8: Index k (fraction of unstable directions) of the most probable saddles as a function of their energy, for the p-model 
(analytical result) and for two models of glass formers (numerical results). In agreement with the mean-field scenario, the most 
probable saddles seem to become stable below a certain energy threshold e t h- For the p-spin model, k oc (e — e t h) i ^ 2 ■ Courtesy 
of A. Cavagna and I. Giardina. 



B. Worrying items 

The most striking feature that appears to be in blunt disagreement with the MCT predictions is that activated 
events are in fact already present well above Td'. 

• the average energy of the inherent structures seem to dip at an onset temperature To > Td jl76l Il77| ; 

• numerical work clearly reveal activation between 'meta-basins' that dominate the long time dynamics above the 
estimated MCT temperature [l78L fl79l |: 

• Vogel-Fulcher fits and Adam-Gibbs correlatio ns w ith the configurational entropy have been reported above Td 
both for experimental and numerical data [20L 1 1 3 ll ] : the power-law increase r Q oc (T — T^) -7 predicted by MCT 
holds at best over three decades in r Q and even this has been questioned; 

• The shape of the relaxation function C(t) in the a region, or of the a peak in dielectric spectra, are suspiciously 
similar above and below Td (see e.g. |l7Q|). 

All these elements suggest that the cross-over at T* ~ Td - if it indeed exists - is not very sharp, whereas the RFOT- 
MCT scenario suggests, at least naively, that the relaxation channels in the a-regime should be completely different: 
unstable modes at high temperatures crossing over to activation at lower temperatures. 

A related worry is the fact that all estimates (numerical and experimental) of the dynamic correlation length ^ 
are remarkably smooth around Td- This is p uzzli ng because the simplest (and probably too naive) interpretation of 
the analytical results of Franz & Montanari [105| suggests that £d should grow as id oc (T — Td)~ v until the MCT 
relaxation time becomes larger than the activated relaxation time related to the mosaic length scale. This should 
occur when: 

id ~ (64) 

where z is the MCT dynamical exponent, given by z = in mean- field (see Eq. (|14p ). But because of the exponential 
dependence of the activated time scale, this crossing is expected to occur whilst i* <C id, suggesting that the dynamic 
correlation length ^ may in fact be non-monotonous with temperature, with a (curious) drop around Td- Of course, 
£d and i* may not necessarily be the same length at all, but it still a challenge to explain how the two regimes patch 
together seamlessly. 

Let us finally note that the value of Td calculated within the MCT approxima tion is usually substantially higher 
than the value needed to fit MCT predictions to experimental or numerical data |l8dj . This could be considered as 
worrying for the MCT approach altogether. However, as shown in [66|, MCT should in fact be understood as a kind of 
Landau theory for the glass transition, that predicts generic scaling form for the relaxation function. In particular, all 
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the predictions quoted in section [II B 31 above are valid beyond the MCT approximation, which might indeed be quite 
imprecise as far as the value of Td is concerned. We believe that a more accurate prediction of T d should be obtained 
using static replica theory instead (see [3] and Appendix A) . The only case where this has been done is hard spheres, 
where the critical dynamic volume fraction (\>d was computed to be <j>d = 0.5159 using the MC T ap proximation and 
4>d ~ 0.62 using replicated HNC, indeed closer to the numerically determined value <j>d = 0.592 [l8l| . 



C. Ways out of the conundrum: a Ginzburg criterion for MCT? 

The presence of activated events above Td are in fact expected for systems of finite size and/or finite dimensions. 
Take for example a mean-field model such as the Random Orthogonal Model (ROM) for which unstable saddles and 
stable minima are strictly "de- mixed" for infinite system size N — > oo. For finite size N, there is always a probability 
(exponentially small in N) to find the system stuck in a minimum above the (sample dependent) transition temperature 
Td- Even if these events are rare, they clearly dominate the long-time equilibrium relaxation C(t). These rare events 
are even more important when one considers the disorder averaged correlation function C(t), because of the order 
A -1 / 2 sample to sample fluctuations of the critical temperature T d : at a given temperature, some samples are 'hot' 
and dominated by the MCT unstable mode scenario, whereas others are 'cold' and dominated by activation over finite 
barriers. The physics of these sample-to-sample fluctuations has been studied in great details for the ROM in [96| 
recently. Numerical simulations of the dynamics of ROM clearly reveal long time relaxation tails above Td that are 
not accounted for by MCT. 

Similar fluctuation effects are also expected in finite dimensions and may affect the MCT scenario in a profound 
way. One can classify spatial fluctuations in three different categories: 

• the usu al cr itical fluctuations that renormalize the MCT exponents below the upper critical dimension d c = 8 

El SGI!; 



• the 'nucleation' fluctuations that lead to activated dynamics; 

• the finite dimensional counterpart of the sample to sample fluctuations found in the ROM, and expected in any 
disordered system. 

Let us focus on the last two. In order to understand their relevance, we first derive a Ginzburg-Harris criterion for the 
MCT transition in finite dimensions. Neglecting for a while activated processes, we recall that the liquid above Td is 
formed by dynamically correlated regions of size i d - This length measures the typical extension of the unstable modes 
and diverges as id oz (T — X^) -1 / 4 . The free-energy of a correlated region, computed within a time scale less than 
T a , fluctuates from region to region. Within the simple REM (parabolic) description of the configurational entropy 
used throughout this paper, these fluctuations are of order Sf(£ d ) ~ T g JEC~ v l d d/2 , which can be interpreted as local 
temperature fluctuations, of order: 

ST(td) = ^ J_ (65) 



dT 



T g JAC p £« 



Assuming mean-field is correct, id oc (T — Td) VmS with v m j = 1/4 one finds that the temperature fluctuations within 
a region of size id are of order: 

6T(£ d )~jl^!=(T-T d ) d / s . (66) 

Clearly, 5T(id) should be much smaller than T — Td itself, otherwise the nature of the slowi ng d own would be totally 
changed. Whenever d > 8, this criterion is always satisfied in the critical region T — > Td [230]. In d < 8, however, 
this is only valid outside a Ginzburg-Harris region defined by: 



T d \T g / ACp-" 

For MCT predictions to hold, on the other hand, one must have T — Td <C Td- Both criteria can be simultaneously 
satisfied only if AC P 3> 1, i.e. for very fragile systems. The dynamic length corresponding to the boundary of the 
Ginzburg-Harris region reads: 

2 

Igh ~ l^r) ACr d (68) 
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which for d = 3, Td/T g ~ 1.3 and AC P ~ 20 (corresponding to very fragile systems, see Fig. 3-b) is still very modest: 
Iqh ~ 1-65! The corresponding relaxation time 

log 10 ^ = 4 7 log 10 ^ Gff w2 (69) 

TO 

for 7 = 2.5. This could explain why MCT predictions at best only explain the first 3 decades of increase of the 
relaxation time, before crossing over to a new regime. 

Within the Ginzburg-Harris region, fluctuations are important and the exponent governing the divergence of Id 
must be renormalised from v m f = 1/4 to a val ue v > 2/d, such that the above criterion is asymptotically satisfied 
(y > 2/d), or marginally satisfied (y — 2/d), see [1161 ] for a proof that this inequality must be satisfied generically for 
disordered systems. 

Let us now consider the fact that there are other types of fluctuations as well. The usual critical fluctuations should 
only change the value of the exponents close enough to Td- In this case the above arguments above still hold but with 
a value of v possibly different from the mean field exponent v m t = 1/4 if the standard Ginzburg region turns out to 
be wider that the above Ginzburg-Harris region, which is unlikely unless AC P is really large. Fluctuations related to 
activated events, on the other hand, are in competition with the ones above, induced by the self-consistent disorder. 
These are the only ones remaining in d > 8 and they also blur the MCT transition. Understanding the interplay 
between these two types of fluctuations is an open problem. We sorely lack deeper analytical investigations of the 
MCT transition in finite dimensions, and we can only venture to formulate some uncontrolled conjectures. 

A possible scenario is that v is precisely equal to the lower bound 2/d; which means that the probability that 
T — 5T(£d) < Td is always of order unity, such that Td cannot be sharply defined. In this case the MCT transition 
is blurred in the following way: a finite fraction 1 — <j> of regions of size Id are governed by the MCT relaxation 
mechanism and "fast" , while the complementary fraction <fi are effectively below the dynamic transition. These latter 
regions have lost all unstable relaxation modes and should be governed by the slow mosaic relaxation mechanism 
with £*(Td) ~ 1. If this is the case, the short time part of the relaxation function C{t) > <fiq* is dominated by the 
faster, MCT channel whereas the long time part, corresponding to C(t) < (f>q* is dominated by activation processes. 
In fact, the asymptotic relaxation regime based on the droplet argument above [see Eq. ([53"]) ] should still hold above 
Td, because of the unavoidable presence of rare, but deeply stable glassites up to the onset temperature To. This 
naturally explains why activated effects are indeed observed above Td and why a peaks look rather similar above and 
below Td- 

Conversely, there should still be occasional unstable soft modes below Td, that dominate the short time behaviour 
of the correlation function. It would be very important to formulate these arguments more precisely, even at a 
phenomenological level. In any case, the following ideas seem to us crucial to understand the MCT/RFOT crossover 
and avoid contradictions with empirical data: 

• Fluctuation effects dominate the physics around Td and lead to a coexistence of MCT and activated dynamics, 
which are hard to disentangle on intermediate time scales. This does not happen in mean-field situations, nor 
in the Kac limit, and requires analytical tools unavailable at this stage. 

• In the temperature-time plane, there should only be a sliver region that is well described by MCT predictions 
(see Fig. ©. 

• The crossover between the MCT relaxation mechanism and the activated mechanism in practice takes place 
whilst all correlation lengths £d, and £* are still very small and close to each other such that no 'kink' or non 
monotonicity of £d(T), that could indeed occur in mean-field, will ever be detected experimentally. 

The above arguments also suggest an interesting lead to study the MCT transition in finite dimensions: the important 
parameter that controls the width of the Ginzburg-Harris region appears to the specific heat jump AC P ; at least one 
kind of fluctuations can be tuned down for models systems such that AC P — ► oo. This may allow one to understand 
the MCT cross-over in more details. 



VI. COMPARISON WITH OTHER THEORETICAL APPROACHES 

In order to assess the achievements of RFOT and understand its strengths and weaknesses, it is useful to compare 
the basic premises of the theory, as well as its most salient predictions, to those of other approaches. There is in fact 
a bevy of different theoretical pictures for the glass transition. We certainly do not want to review all of them but will 
pick of few proposals that seem relevant to us, either because they rep resent viable alternatives to RFOT, or because 
they have been actively discussed in the recent years (or both !). [23 1| 
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FIG. 9: Cartoon of the conjectured MCT-RFOT crossover in the time-temperature plane for d < 8. We parameterise time 
between and 1 as 1 — C(t) (x-axis). Below the hypothetical Tk, C(t — * oo) = q* > 0, and the red region is inaccessible. The 
yellow region corresponds to activated dominated dynamics, that extends up to the onset temperature To > Td if t is large 
enough. The blue MCT "sliver" region is confined to the immediate vicinity of Td and for short times, mostly in the j3 regime. 
The green region at high temperatures corresponds to free flow, when barriers are irrelevant because typical saddle points are 
highly unstable. 



A. Elastic models 

Dyre's "shoving model" [l46l Il47| postulates that elementary events that allow supercooled liquids to flow are 
nucleation of "voids" that allow particles to move around and unjam. The energy cost A of such voids comes from the 
elastic deformation of the surroundings during the short instant during which the void is created. This deformation 
only involves shear, and Dyre argues that the corresponding energy barrier is: 

A(T) cx G^TjJZ* 3 (70) 

where R* is related to the size of the critical void, which is assumed to be temperature independent, and Goo the high 
frequency shear modulus, relevant for these rare, but supposedly quick, activated events. 

The non-Arrhenius behaviour of r Q is therefore entirely ascribed to G^; the crossover temperature T* is associated 
to a noticeable stiffening of the liquid at high frequencies. Since the thermal fluctuations of the particles around 
their amorphous equilibrium position is also dominated by shear modes, the plateau value u 2 of the mean-squared 
displacement are also given in term of G^ as: 

2T 

u 2 p = ( U 2 (r «t«r Q )) » Q-^j- ( 71 ) 

Interestingly, the two equations above relate the short time displacements of the particles to the slow relaxation time 
of the system. Eliminating A(T) leads to: 

^~=A- 2 , (72) 

TO 

where a is the average inter-atomic distance and A is a tem perature independent numerical constant. This relation 
was actually first proposed by Hall & Wolynes in 1987 jl86| . but within a framework quite different from RFOT - in 
fact, this relation is very unnatural within the context of Wolynes' "standard" version of RFOT, see below for more. 
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Now, there are a number of very good points about Dyre's theory. First it is extremely simple and the assumptions 
are straightforward and transparent. Second, the model fares quite well at accounting for the non-Arrhenius depen- 
dence of r a with temperature. Plotting lnr Q as a function of Goo/7 1 l eads to much better (but not perfect) straight 
lines than when plotted against 1/T, as indeed predicted by Eq. (|TD|) [l44l 1 1451 . Il47| . Third, the correlation between 
In 77, and short time motion, Eq. (|72[) , is well supported by recent experiments [191 ]. Actually, recent simulations 
[l87l | have shown that this relationship holds even locally in model supercooled liquids: 'soft' regions where Up is large 
indeed relax faster than 'hard' regions. 

In this purely local picture, where the super- Arrhenius behaviour is totally determined by Goo, the Adam-Gibbs 
correlations is deemed not fundamental but merely fortuitous. Despite this fact, the model, if correct, has at least 
to be compatible with these Adam-Gibbs correlations: between energy barriers and excess entropy on the one hand, 
between jump of specific heat and fragility on the other hand. A way out is suggested below. Anothe r cor relation 
that should be explained is why more fragile materials seem to have a lower stretching exponent j3 0, Il5dj | , and a 
stronger violation of SER [4(J, suggesting more fluctuations. Of course, one can always argue that these correlations 
are, again, nonexistent or artefacts. 

In summary, the correlation between the effective energy barrier and the high-frequency shear modulus, that explains 
at least part of the super- Arrhenius behaviour, is certainly striking and suggestive. The most serious issue with this 
approach is, however, related to its core assumption: that the glassy slowing down is a purely local phenomenon, 
and that the size R* of the critical void is temperature independent. This is at odds with our current understanding 
of dynamical heterogeneities, and with all recent theo retical and numerical results on the growth of a characteristic 
length as temperature is decrease d (see e.g . \I§L Ii~74l0 . Experimental results on dynamical correlations also clearly 
favour a growing length scale [iol. Il38l 1 139L ll42Ul43l |. In particular, the very recent measurements of the non-linear 
dielectric constant of glycerol unambiguously reveal tha t coherent amorphous order indeed propagates over larger 
and larger distances as the temperature is reduced jl4l| . Although one could still argu e that the se facts are mere 
consequences and not primary causes of the slowing down, recent numerical simulations [l2dl . Il35| seem to rule out 
that the growth of the effective barrier is due to a purely local process. In particular, as we mentioned in the discussion 
of finite size effects above (see section HV Ap , the relaxatio n tim e r(i?) of a cavity of size R grows with R up to the 
thermodynamic mosaic length £* and saturates for R > £* [l35| . If the explanation of the slowing down was local in 
space and only due to the growth of Goo(T), one would expe ct no change in t(R) as soon as R > R* , where R* is 
temperature independent, contrary to numerical observation. [2321 ] Note that since Goo(T) is expected to saturate at 
low enough temperatures, the energy barrier should stop increasing if R* is indeed temperature independent. In the 
Dyre model, r Q should revert to a purely Arrhenius behaviour at low temperature, as indeed expected if no growing 
order of any kind is present in the system. Finally, a p ossib le problematic point of this model is the value of R*, found 
to be in the range 3 — 4 Angstroms using the data of 145]. It would be very important to rationalise this number in 
terms of a realistic microscopic rearrangement process. 

Dyre's model and RFOT might in fact be more akin than may appear at first sight. Even if one accepts that an 
important contribution to the super- Arrhenius behaviour is indeed the anomalous growth of Goo(T) when T is less 
than the Goldstein temperature T* , then we should explain the microscopic mechanisms for such a sudden change of 
behaviour of the system around T*. Here, the Goldstein argument of a crossover between unstable saddles and stable 
minima seems unavoidable, because any growth of Goo is nearly tautologically related to the stability increase of the 
local configurations. Therefore, an MCT-like mechanism should be at play, at least to understand the behaviour of 
Goo{T). Within an 1-RSB like scenario, the increase of G 00 (T 1 ) is associated to a change of dominant metastablc 
states: as T is decreased, the system probes deeper and more stable minima of the energy landscape. Interestingly, 
as noted in section IIII El above (see also Fig. [5]) , the rapid change of Goo (T) is in fact associated to an important 
change in the vibrational entropy. This could provide a natural explanation of the Adam-Gibbs correlation, without 
any direct link between the configur ational en tropy and energy barriers, but only a co-variation of these quantities, 



But if this 1-RSB mechanism is indeed the explanation for the increase of Goo(T'), the mosaic argument should be 
valid too, re-introducing the increase of the mosaic length I* as a source of super- Arrhenius behaviour. A possibility 
is that in the experimental temperature range, the energy barrier A(T) = increases both because Ao oc Goo 

increases, as in Dyre's shoving model, and because £* increases, although in a very modest fashion. This would allow 
one to account for the clear curvature in the plots of lnr a as a function of Gqo/T (see e.g. [l45| . Figs. 5 and 6), 
and reconcile Dyre's mechanism with the existence of a growing length scale. On the other hand, all the quantitative 
estimates based on Wolynes's choice of RFOT parameters (see section HV Bp would need a complete overhaul. Another 
possibility is that the barrier crossing leading to relaxation on the scale £* involve smaller and smaller displacements 
per particle, 8, upon lowering the temperature. This would lead to an effective barrier that is of the order Goo^ 3 ^*^, 
where 6 3 decreases approaching Tr- and l*^ increases, providing an interpretation of the correlation between A(T) 
and Goo, and leaving room for £* to grow. In any case, a convincing reconciliation between the shoving model and 
RFOT is needed for both theories; this means in particular elucidating the relation between £* and R* (if any). This 
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seems to us a very important point to clear up in future research. 

Finally, coming back to the issue of a growing length scale, a clear-cut distinction between the shoving model with 
a fixed R* and the mosaic picture with a growing £*(T) is the behaviour of the cross-over stress a* beyond which 
shear-thinning is expected. In the mosaic picture, a* should scale as ^JTG^TJ^T). In the shoving picture, on the 
other hand, a* is simply given by Goo- This should lead to testable differences: a relatively fast growth of a* when 
T decreases for the shoving model and a roughly constant or even a mild decrease for RFOT. 



B. Frustration Limited Domains 

This class of models has been advocated by Kivelson, Tarjus & c oworkers s i nce t he mid-nineties [HI, HH, [23[ , building 
upon older ideas dues to Frank, Sadoc et al., Nelson and others jl89l Il90l Il9l| . The basic tenets of this approach 
have been very clearly articulated in the review which we are satisfied to reproduce almost identically here: 

1) A liquid is characterised by a locally preferred structure (LPS) which is different than that of the crystalline 
phases. 

2) Because of geometric frustration, the LPS characteristic of a given liquid cannot tile the whole space. 

3) It is possible to construct an abstract reference system in which the effect of frustration is turned off. 
Frustration naturally leads to collective behaviour on a mesoscopic scale (■■■)■ The collective property comes from 

the growth of some ordered phase (in which the liquid LPS spreads in space) induced by the proximity of an (avoided) 
critical point. The limitation of the scale over such a growth can take place results from frustration that aborts the 
phase transition and leads to "avoided criticality" . 

There are some analogies, but also deep differences, between the ideas of Frustration Limited Domains (FLD) and 
RFOT. As in RFOT, there is a crossover temperature T* ~ T4 where some kind of local thermodynamical order set 
in: amorphous TAP states in one case, PLS in the other. This gives rise to high energy barriers because collective 
motions are needed for any rearrangement to take place. As in RFOT, the size £* of the glassites are limited for any 
given temperature T < T*, and grows when T decreases. As in RFOT, the energy barrier is postulated to grow as a 
power of £*. 

However, the size £* is limited by configurational entropy in the case of RFOT, and by long range elastic energy 
that frustrates the local order in the case of FLD. More precisely, the proponents of FLD argue that £* grows 
as (1 — T ' )T*Y I 'K 1 / 2 where K is an a-dimensional parameter measuring the strength of frustration, and v is the 
exponent governing the growth of the correlation length of the unfrustrated transition (K = 0). As expected, £* 
diverges if frustration is absent (K — ► 0). The energy barrier for the relaxation of the glassites is assumed to be given 
by A = AT£* 2 , where T(T) is the surface tension between two different LPS. Since the (avoided) phase transition 
towards LPS takes place at T*, it is reasonable to assume that T(T) itself vanishes at T* as T oc T* 1 " 2,7 (T* - T) 2v . 
This leads to an energy barrier that behaves as [H, H2, HH : 



A(T<T*) = A > + — (^1- — J , (73) 

where we have added a high temperature energy barrier that governs the frequency of elementary moves. This has 
the shape alluded to earlier, see Eq. (|63|) and was used to fit the data shown in Fig. 2 with v in the range [1/2, 3/4] 
[22, The value v = 2/3 in fact corresponds to an Ising like transition of the unfrustrated model in d = 3. The 
model was extended to account for fluctuations of £* and predict the broadened dielectric spectrum as a function of 
temperature and frequency, with remarkable success [l2l |. 

However, we note that below a certain temperature, sufficiently far from T* , both £* and T should saturate and the 
barrier A(T) should become constant. As within the shoving model, FLD predicts reversion to a purely Arrhenius 
behaviour at low temperature, whereas RFOT suggests that A(T) continues to grow as the configurational entropy 
goes to zero. 

Recently, this scenario was carefully tested 192] by simulating a mono-atomic liquid on the hyperbolic plane, where 
crystallisation is avoided because the non-zero curvature frustrates the crystalline order. As a consequence, this setting 
provides, mutatis mutandis, a benchmark to study the FLD theory. The results that have been found are in overall 
agreement with the above s calin g arguments and assumptions. Signatures of LPS in glass-forming model systems 
have also been reported, see [l93j . 

There are two main issues that we would like to address concerning FLD theory. The first one is its compatibility 



with the Adam-Gibbs like kinetic/thermodynamic correlation. Kivelson and Tarjus [194j have argued that FLD could 
cope reasonably well with the correlation between excess entropy and relaxation time. What is unclear to us is whether 
the correlation between fragility and specific heat is compatible with FLD. Simple - but arguably not very reliable 
- arguments suggest that there might indeed be a problem. The main difficulty is to obtain a good estimate of S xs 
within FLD. 
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Kivelson and Tarjus argued in [194j that the proximity of the unfrustrated phase transition leads to a singular 
contribution to the excess entropy given by S xs w S(T*) - A(l-T/T*y u \ where A is a numerical constant 
independent of K, and hence of the fragility. [233| | Although this expression indeed predicts a decreasing S xs below 
T* , it implies an incorrect evolution with fragility, since it leads to a specific heat jump at T g totally independent of 
fragility! 

Another physically transparent (but still moot) estimation the configurational entropy is the following. Let Sq be 
the total zero temperature entropy of the unfrustrated model (i.e. a measure of the degenerescence of the LPS). 
Since each glassite of size £* can independently be in e So different LPS, the configurational entropy per particle at 
temperature T is given by: 

E(T) = || oc S K 3/2 (l - T/T*)-^. (74) 
Comparing with Eq. (|73|) . this leads to an Adam-Gibbs relation for a given material (i.e. a given K, So)' 

, r a A> A .T*st^K , . 

with an exponent —4/3 instead of the usual Adam-Gibbs exponent — 1. These values are probably hard to distinguish 
experimentally. 

Now, taking the derivative of the above relation with respect to temperature allows one to compute the fragility to 
as: 

(76) 



m 3£(T g ) T* - T, 



'i 



where we neglect the contribution of A> at T g for fragile materials, which turn out to be such that K — > 22]. The 
above direct relation between fragility and T* — T g could be tested empirically. This would allow one to measure v 
directly. Within the same approximation, the very definition of T g leads to: 

AT*(1-T g /T*)^ T* ~T g (m Q K\ 1/Au 



T g K T* 



showing indeed that to — > oo when K — > 0. From this relation, we also find that the size of the domains at T g is 
£*{T g ) cx K- 1 / 4 and (weakly) grows as a function of fragility, whereas S(T ff ) oc SoK 3 ^ 4 . The last information we need 
is AC P computed from the explicit form of S(T): 

ACp * 3y (l-T°/T*)^ * ^ o(1 " T s/ T *) 3V ' 1 « S.K 1 ^, (78) 
which, combined with the above relation finally leads to: 

i_ 

m — m ( AC p \ 1-3 



oc =P ■ (79) 
m \ So J 

If we assume that So does not vary much between different materials, while the change of K is responsible for 
the change of fragility, we find that m and AC P should be inversely correlated whenever v > 1/3. For example 
(to — TOo)/m oc ACp 1 for the favoured value v = 2/3. Another way to see this is that AC P goes to zero when K — > 
when ^ > 1/3. The previous estimation of the configurational entropy does not take into account the entropy due 
to the wandering of domains between LPS. If So is not a constant but scales with £* then E(T) would also have a 
different scaling with £*. Assuming S(T) oc l/£* 3 ~ d J , where df is a certain fractal dimension, one can find the correct 
correlation between specific heat jump and fragility for some range of df. However, for these very same values of df 
the Adam-Gibbs relation would be completely altered in a way detectable in experiments. 

As a conclusion, it seems that obtaining the correct evolution with fragility of the configurational entropy and the 
specific heat jump is a very non trivial constraint for FLD theories. Because of the absence of an analytically reliable 
way to estimate the configurational entropy it would be worth addressing this question using simulations, in particular 
on the hyperbolic plane where FLD ideas have been shown to be valid [192J. 

Another quite interesting aspect of FLD theory is that any reasonable analytical approximation used to deal with 
explicit models for which the FLD scenario is expected to hold (such as, for example, Brazovskii's frustrated Coulomb 
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model) in fact all lead to RFOT-MCT equations @, [H, HH HH ■ The crucial assumption that the number of LPS of 
the system is small, and not exponential in N, could therefore be unwarranted. An exponential number of amorphous 
metastable states emerge instead, induced by frustration. This extensive configuration entropy is in fact at the root 
of the positive correlation between fragility and jump of specific heat, captured by RFOT. 

Finally, the crossover stress in the FLD theory should be given by a* oc \jGoaT /I* ~ K 1 / 4 (1 — T/T*y/ 2 , which now 
increases as T is decreased. When T = T g , a*{T g )/y%G^ should scale as K 3 / 8 for very fragile materials. Actually, 
the FLD behaviour of a* when T is lowered is intermediate between the shoving model, which corresponds to an 
increase proportional to Goo, and RFOT, which corresponds to a variation slower than y/Uoo. Non- linear rheology 
data for a range of fragile materials would certainly be very useful to distinguish between all these pictures. 

C. Others 

Another active line of research is devoted to Kinetically Constrained Models (KCM, for reviews see e.g. [36l . Il95| |). 
The idea behind this family of models is that one might be able to find a coarse-grained effective description of the 
dynamics of glass forming liquids where the complicated microscopic interactions that slow down the dynamics are 
replaced by simple "kinetic" constraints. For example, one could replace the initial problem by a lattice model, where 
each cell contains a certain number of particles, enough for the coarse-grained description to be valid. Each cell is then 
"active" if some motion can take place, or "inactive" if jammed. The density of active cells is assumed to be given 
by a simple exponential: p = cxp(— J/T), where J is an activation energy. The number of active cells is supposed to 
be strictly conserved: activity cannot spontaneously appear nor disappear, but only hop to nearest neighbours. This 
is where kinetic constraints come into play. A simple rule is that the hop is only allowed if the target cell is itself 
surrounded by a sufficiently large number of active cells. Several other rules have been considered, leading to different 
models, some exhibiting very int eres ting properties, such as the existence of a genuine dynamical transition at some 
critical density of active cells p c jl96l |. 

But the truly interesting outcome of these models is that kinetic constraints lead to cooperative dynamics: mobile 
regions must cluster together in order for the system to evolve. Without constraints, the relaxation time r a would be 
simply proportional to p~ x in the dilute regime p — > 0, leading to an Arrhenius slow-down. By imposing cooperativity, 
kinetic constraints can deeply affect this behaviour, and give rise to super- Arrhen ius b ehaviour, for example r Q ~ 
p~ J / T as T — > 0, which reproduces Bassler's 1/T 2 law mentioned in se ction II Bl [ 36l. Il95l | . KCM also offer interesting 
ideas to rationalise dynamical heterogeneities and SER violations |l48| . 

However, the basic assumption here is that kinetic constraints can lead to non-trivial and interesting effects, even 
if the thermodynamics is completely trivial. By construction, KCM has a hard time explaining correlations of the 
Adam-Gibbs type. In fact, the large specific heat jump AC ? at T g for fragile glasses (several ks per particle), cannot 
be accounted by the freezing of dilute mobility defects |l97| . The thermodynamics of each cell in the coarse-gr aine d 
model must be contributing to AC P in order to cure the problem of having too small a specific heat jump at T g 198]. 
On the other hand, the correlation between kinetic and thermodynamics, and its evolution with fragility, then requires 
that these correlations are encoded in the thermodynamics of the cells. But this appears to be self-contradictory, 
because by assumption nothing at all happens within inactive cells, while these represent the overwhelming majority 
of cells as T approaches T g - so nothing big should freeze at T g since it is already frozen above T g \ 

In any case, the crucial issue is to work out how a system of dense, interacting particles can really be mapped onto 
a KCM at sufficiently large length scales. The mechanism by which the system becomes locally rigid [234| in the 
first place must be accounted for. Most people, including the advocates of KCM, now seem to agree that this local 
rigidity sets in around a rather well defined temperature T* (see e.g. [13])) this is an essential ingredient of all theories 
reviewed so far (RFOT, elastic models, FLD and KCM). In order to account for such an abrupt change in the local 
rigidity and the appearance of constrained dynamics below T* , an MCT type of scenario seems to us difficult to avoid 
- with stable, amorphous and/or locally preferred structures, becoming thermodynamically predominant below T*. 

Once this local rigidity is installed, the scale beyond which an effective, simplified description is possible should 
correspond to the size of an elementary cooperative volume, i.e. the number of particles involved in an elementary 
activated event. This would be £*(T) in the framework of RFOT or FL D, an d R* in the shoving model. However, 
this issue has only been very cursorily discussed in the KCM literature jl56l |. although crucial to understand what 
regime of length and time scales KCM are supposed to describe. As mentioned in section llV CI one expects that some 
facilitation mechanism should be present also within RFOT on length scales larger than I* : when an activated event 
takes place within a droplet of size £* , the boundary conditions of the nearby droplet changes, which can trigger a 
second activated event, possibly inducing an "avalanche" process that extends over the dynamic correlation length 
scale £d > £*. However, we fail to se e why activity shou l d be an almost conserved quantity, as assumed in the series of 
papers of Garrahan & Chandler [H EM EH EM EH, GUI ■ Actually, we believe that a more likely situation is that 
activity becomes less and less conserved upon lowering the temperature, as it happens in granular systems close to 
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the jamming transition [157l . ll58j |. Spontaneous, thermally activated events occurring within immobile regions should 
be more and more important as the system slows down. Clearly, this is a very important issue that should hopefully 
be settled by numerical simulations. 

Finally, we want to mention the work of M. Moore & collaborators [92. lolil Il99l l200j ] . where the working hypothesis 
is that all models of glasses that are 1-RSB in mean field can in fact be mapped onto a spin-glass in a magnetic 
field in finite dimensions. This is indeed supported by some analytical and numerical calculations. Although some 
aspects of the physics of spin glasses look very similar to that of glasses, there are also striking differences, such as the 
continuous vs. discontinuous nature of the order parameter at the transition. Nevertheless, some efforts to reproduce 
the phenomenology of glasses within that framework have been reported [200j | . 

If Moore's conclusion turned out to be valid, this would be quite a fatal blow for RFOT, except if the time 
and length scales below which the mean-field 1-RSB physics survives are large enough to cover the experimentally 
accessible wi ndow . However, some of the calculations bolstering this strong claim are quite complex and rely on subtle 
assumptions |l99j . which in our opinion are not yet firmly established. Moore's conclusion is also at odds with the 
success of 1-RSB calculations for realistic models of liquids and hard-spheres reviewed in Furthermore, recent 
numerical simulations of some finite dimensional models clearly de mon strate that 1-RSB effects do seem to survive 
away from mean-field, at least over a large time window [U [97], 1 2 ll | . This whole issue is related to the fate of 
MCT in finite dimensions, that we have commented on in section IV CI We lack firmly grounded theoretical work on 
the interplay between 1-RSB and finite dimensional fluctuations, without which RFOT will remain a beautiful, but 
fragile, phenomenological theory. 

VII. DISCUSSION AND CONCLUSIONS 

A. Successes and Difficulties of RFOT 

Let us summarise the strong and the weak aspects of RFOT that are strewn in various parts of this paper. First, 
some undeniable successes: 

• From a theoretical point of view, RFOT is a natural finite dimensional interpretation of the mean-field theory of 
generic complex systems with many minima. Two characteristic temperatures appear naturally: (i) a dynamical 
temperature where the phase space gets fragmented into different minima, which in mean-field corresponds to 
a Mode-Coupling non-ergodic transition, and in finite dimensions to the temperature where some high-frequency 
rigidity appears and leads to activated dynamics; (ii) a static temperature Tk where the configurational entropy 
vanishes and where the system, if equilibrated, would enter an ideal glass state. These two temperatures echo 
the well-known Goldstein and Kauzmann temperatures that describe the phenomenology of supercooled liquids. 
Item (i) appears to be a crucial ingredient of all viable theories to date (RFOT, FLD, elastic models, and even 
KCM): glasses lose their rigidity as temperature increases because at some point local stability is lost. 

• Any reasonable analytic approximation that describes the thermodynamics or the dynamics of realistic models 
of supercooled liquids, using e.g. density functional techniques, replica theory, the Bethe-Peierls approximation, 
the projection operator formalism or self-consistent re-summation schemes, all lead to a 1-step replica symmetry 
broken low temperature phase (characteristic of systems with an extensive configurational entropy) and to a 
non-ergodic transition of the type predicted by Mode-Coupling Theory. Mode-Coupling Theory can actually be 
derived within a Landau approach of discontinuous glass transitions, as a generic expansion where the analogue 
of the 'order parameter' is the (small) difference between the correlation function C(t) and its plateau value 
<f{Td) > 0. From all this, it is hard to see how a truly first principle statistical mechanics approach to the glass 
problem can lead to a very different theory. It is always possible that strong fluctuations or non perturbative 
effects completely change the picture in three dimensions, as is the case for example of the Kosterliz-Thouless 
transition in two dimensions. However, because of the firm mean-field foundation of RFOT, one would have to 
come out with compelling physical arguments for such a scenario. 

• RFOT suggests that the liquid between Tk and Td can be thought of as a mosaic of glass nodules or "glassites" 
with a spatial extension t*(T) limited by the configurational entropy. Regions of size smaller than 1*{T) are 
glasses: they are effectively below Tk and cannot relax, even on very long time scales. Regions of size greater than 
I* are liquid in the sense that they explore with time an exponentially large number of unrelated configurations, 
and all correlation functions go to zero. The relaxation time of the whole liquid is the relaxation time of glassites 
of size t* . The crucial assumption, that thermodynamics alone fixes the value of £*, is well borne out by static 
and dynamic simulations of a liquid inside a cavity. The existence of a thermodynamic "point-to-set" correlation 
length t* that diverges as T — > T£ is supported by exact analytical calculation in the Kac limit of weak, long 
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range interactions |105|. This scenario is also compatible with recent experiments that measure the non-linear 
(third harmonic) dielectric constant X3( w ) °f a standard molecular glass, glycerol. If some transient amorphous 
order sets in over a length scale I* , X3(^) is expected to grow as a power of £* when wt„ ~ 1 but to remain 
small for u> — > 0, since no static transition takes place, at variance with sp in-gl asses where X3( w — 0) diverges 
at the transition. Precisely such a behaviour was seen in the Saclay group [141| . 




• Within RFOT, the large excess entropy of supercooled liquids at T g (several ks per molecules) is attributed to 
the existence of an exponential number of unrelated stable amorphous configurations. The essential role played 
by the configurational entropy allows one to account very naturally, in a way not too sensitive to additional 
assumptions, for a series of empirical correlations between thermodynamics and dynamics: Adam-Gibbs relation 
between the relaxation time and the configurational entropy, correlation between the fragility and the jump of 
specific heat AG P , and between the stretching of the relaxation function and AC P . At least some of these 
correlations are not addressed, or difficult to interpret within alternative descriptions, such as Dyre's shoving 
model, the Frustration Limited Domain theory, or Kinetically Constrained Models. Although one can of course 
take the view that all these correlations are entirely accident al, w e do believe that they impose strong constraints 
on the theory of glasses and have to be taken into account [235| . 



There is also a worrying list of loose ends, caveats and difficulties: 

• The existence of the mosaic length £* heavily relies on the hazy concept of "surface tension" or mismatch energy 
between different amorphous states. It is not easy to give a pre cise oper ational meaning to this idea, that 
would allow one to compute or measure this quantity (see however |l06l . Il07| for an attempt) . Correspondingly, 
there is no consensus, neither on the value of exponent 9 determining the size dependence of such a mismatch 
energy, nor on the dimensional prefactor To and its temperature dependence. The situation is even worse for 
energy barriers. We have no idea whatsoever about the nature of the activated events allowing the system 
to jump between on amorphous state to another on scale £* - are these made of compact droplets, fractal 
objects, excitation chains, void nucleation? The exponent ip and prefactor Ao fixing the height of the barriers 
as A = Aq£*^ are even less understood. Wolynes' prescription corresponds to 9 = ip = d/2 and Ao = To = nT. 
Although based on suggestive arguments, it is by no means obvious and not supported at this stage by any 
first principle calculation. When available, these first principle calculations (analytical or numerical) suggest 
9 = d — 1, but fluctuations effects are probably neglected. We also find difficult to fathom why Ao should 
decrease with temperature. 

• Related to the previous point, shouldn't the energy barrier scale Ao in fact be related to the high frequency 
shear modulus Goo? This quantity is known to have a non trivial temperature dependence that may account 
for an important fraction of the increase of the energy barrier between T* and T g . This is the claim made by 
Dyre and others, in the context of elastic models where the energy barrier is directly proportional to Goo, and 
therefore also related to the short time fluctuations (u 2 ) cx T/Goo of the particles around their average position. 
But if most of the energy barrier increase is due to the behaviour of the shear modulus, there is little space 
left for the essential prediction of RFOT, i.e. that the super- Arrhenius behaviour is due to the increase of the 
mosaic length £* , itself driven by configurational entropy. We have discussed this issue in more depth in section 
IVI Al understanding whether the empirical correlation between the energy barrier A and the high frequency 
shear modulus Goo is fortuitous or fundamental seems to us one of the most perplexing quandary that needs to 
be resolved. 

• Finally, the crossover between a high-temperature MCT regime and a low temperature activated (mosaic) regime, 
supposed to be a strong selling point of RFOT, is very poorly understood even at a phenomenological level. 
There is certainly no clean crossover between a well developed MCT scaling regime above T* and activated 
dynamics below T* , but rather an intricate situation, perhaps resembling the sketch in Fig. [21 where activated 
events appear well above Tj. We believe that this crossover is dominated by fluctuations, due both to the finite 
dimensional effects that affect the MCT regime and to a broad distribution of local energy barriers. 

In order to put the MCT crossover on a less shaky ground, two pieces of evidence are needed: First, a smoking 
gun experiment or a simulation that conclusively shows the existence of a MCT-like mechanism. This could be 
realized by studying finite size effects and equilibration inside small cavities both below and above Td- Second, we 
need a theoretical approach able to lead to precise predictions in the crossover region. This will require some major 
technical breakthrough. Until then, this potentially informative crossover regime will remain somewhat mysterious 
and inconclusive. 
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B. Topics not addressed here, open problems and conjectures 

There is a number of other properties of glasses that RFOT might be able to account for, but that we have chosen 
to leave out, both for lack of space and because we feel that our level of understanding is not satisfactory. These are 
in particular the very low temperature properties, dominated by 'two- level systems', where quantum effects become 
important. Lubchenko & Wolynes argue that some of these low temperature properties are inherited from the mosaic 
structure as the liquid freezes at T g [111 ] . This is a very interesting topic on which we hope to return in the future. The 
relation of the so -called "Boson -peak" with MCT-like soft modes has also spurred a spree of activity in recent years, 
see in particular jll2l |202| . [203] . We have only rapidly touched upon some aging properties when discussing energy 
relaxation after a quench in section IIVD1 see also UJ. Much more has to be done, in particular concerning two- 
time quantities, fluctuation-dissipation theorem, effective temperatures, temperature jumps and temperature cycles. 
One exciting idea is that the notion of "temperature chaos" makes sense in systems without quenched disorder. More 
precisely, one could expect that the dominant state a with in a cavit y with a fixed boundary condition abruptly change 
as temperature is varied, as a result of level crossing (see [2041 .1205]). This could induce some interesting rejuvenation 
effects. 

Although RFOT appears to contain all the ingredients to understand the strong violations of SER in fragile liquids 
(see section HV Cp . we lack a more detailed understanding of the so-called fractional Stokes-Einstein relation that 
approximately holds between T g and T*, i.e. D s r//T oc r„, where D s is the self-diffusion constant and k > an 
exponent that appears to increase with the fragility of the liquid. The argument must be that some activated events 
of size £ < I* contribute to the diffusion constant but not to the stress relaxation, for example by allowing the 
permutation of particles without structural rearrangements, as in a perfect crystal for example. A related poorly 
understood problem is the dynamics of a tracer particle driven by an external force in the disordered, but time 
dependent environment created by the surrounding particles (see [206j ] for recent results in the MCT regime). 

We have not dwelled on the subject of dynamical heterogeneities either. It is clear that while RFOT predicts strong 
spatial heterogeneities and temporal intermittency in the dynamics, it is not clear to us how one should compute the 
dynamical correlation length ^ and compare it to various recent numerical and experimental determinations of these 
dynamical correlations. Whereas cooperative activation events take place over the mosaic length £*, it might be that 
facilitation effects leads to avalanches of activity that spreads over much longer length scales, in which case £d 3> £* ■ 
Cooperativity and correlation do not mean the same thing, although it is by no means easy to define cooperativity 
unambiguously, i.e. in a such way that one devise a precise protocol to measure it (see [1571 ] for a recent attempt in 
granular media). 

What do we need at this stage? First, some major analytical progress on a realistic, finite dimensional system where 
RFOT/MCT is expected to hold, to understand in detail the interplay between spatial fluctuations and activated 
events. We need to understand whether the cartoon sketched in Fig. [51 that delineates the different regions of 
applicability of MCT and RFOT, makes sense or not. The Ginzburg-Harris criterion we discussed above tells us 
that the limit of an infinitely large jump of specific heat should ease some of the problems. A toy-model where non 
perturbative activated effects can be studied thoroughly would be extremely valuable. Some of the ideas sketched in 
the Appendices might turn out to be useful. 

Short of theoretical ideas, one should at least try to find a finite dimensional version of a 1-RSB mean-field model 
that can be convincingly simulated both in the MCT regime and in the activated/mosaic regime. This would allow 
one to test directly some of the building blocks of RFOT. Unfortunately, the only models in that category can only be 
simulated above the MCT temperature T^, and cannot be equilibrated below (at least up to now). Although the study 
of the high temperature side of T e i brings a decent amount of information (in particular on the role of activated events 
in that region) the situation is frustrating. It might be the case that models that are 'hard' enough (in the sense that 
q* = C(to <C t <C r a ) is close to unity) so that 1-RSB effects survive in low dimensions, are also, for the very same 
reason, extremely difficult to equilibrate numerically |96[. Intuitively, these are systems with a "golf-course" energy 
landscape and very narrow canyons. Is this a necessary prerequisite to form a glass? 

One could also hope for a smoking gun experiment of some sort, that either definitely rules out or strongly support 
the RFOT scenario. We have seen, for example, that the role of shear on viscosity is quite different in RFOT and 
in alternative theories, in particular the behaviour of shear thinning as a function of temperature. Energy or volume 
relaxation after a deep quench and other aging/rejuvenation experiments could also be very informative. Any direct 
proof of the existence of a configurational entropy dominated point-to-set length would be decisive. One could surely 
do with more "cavity" numerical simulations and/or real experiments on the statics and dynamics of a liquid trapped 
within walls made of the same frozen liquid. Certainly, the idea that one can force a small enough system of size < £* 
into an ideal glass phase even when T > Tk is worth exploring experimentally. 

Finally, if the relaxation time r a really varies exponentially with the size of the glassites £* , we must come to terms 
with the fact that cooperative regions are doomed to remain small, maybe 10 3 molecules if we are lucky, or £* ~ 10 
at most. Scaling relations are at best guiding lights to understand general trends, but can never be tested accurately. 
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Sub-leading corrections will always soil the asymptotic values of the exponents that toil and sweat might be able 
to produce once quite formidable theoretical challenges unravel. The glass transition problem will probably only be 
solved like in the game of Go, by a slow, patient siege. Because our plea for a smoking gun experiment will probably 
never be granted, we should at least remember Feynman and put down all the facts that disagree with our pet theory 
as well as those that agree with it. 



Acknowledgements 

We want to warmly thank V. Lubchenko and P. G. Wolynes for asking us to write up our thoughts on RFOT, 
and for their patience. We thank L. Berthier, A. Cavagna and G. Tarjus for a careful reading and feedback on the 
manuscript. Finally we also want to thank all our collaborators on these subjects, namely: C. Alba-Simonesco, A. 
Andreanov, L. Berthier, E. Bertin, A. Billoire, R. Candclier, A. Cavagna, L. Cipelletti, L. Cugliandolo, O. Dauchot, D. 
S. Fisher, Y. Fyodorov, T. Grigera, W. Kob, J. Kurchan, F. Ladieu, F. Lechenault, A. Lefevre, D. L'Hote, M. Mezard, 
K. Miyazaki, C. Monthus, D. Reichman, T. Sarlat, G. Tarjus, M. Tarzia, C. Thieberge, C. Toninelli, P. Verrocchio 
and M. Wyart. We have benefited from enlightening conversations over the years with C. Brito, C. Cammarota, M. 
E. Cates, D. Chandler, C. Dasgupta, J. C. Dyre, S. Franz, J. P. Garrahan, P. Goldbart, A. Heuer, R. Jack, J. Langer, 
A. Montanari, M. Moore, S. Nagel, G. Parisi, S. Sastry, P. Sollich, T. Witten, P. G. Wolynes, H. Yoshino and F. 
Zamponi. 

We were partially supported by ANR DYNHET 07-BLAN-0157-01. 



APPENDIX A: ANALYTICAL APPROACHES TO METASTABLE STATES AND CONFIGURATIONAL 

ENTROPIES 



1. Density functional approaches 



As we discussed in section H-B, the solution of some mean- field spin-glasses unveils a very rich and interesting 
scenario that bears some resemblance with the physics of glass forming liquids. Still, there seems to be a long way 
before making a quantitative connection with the behaviour of realistic model of supercooled liquids. The aim of this 
Appendix is show how this goal can achieved, at least to some extent. The specific issues we address are: 

• How can one identify metastable states theoretically? 

• How can one compute their physical properties? 

• How can one check that a given system displays indeed many metastable states (without resorting to a dynamical 
or numerical analysis)? 

For simple systems, like a ferromagnet, the standard procedure consists in identifying the minima of the free energy 
functional. In the case of liquids, this would correspond to the density functional defined as 



n / . r jv 

^({p(x)» = -^log^y / Y[d Xl exp [-0H+ dx M (x) ]T S(x 
" ' J i=l \ U=l 



Xi) - p(x) 
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where H is the Hamiltonian of the system and /i(x) is a chemical energy fixed by requiring that (X)i=i <^(x— Xj)) = p( x )- 
The average in the previous expression is obtained using the Boltzmann weight on the RHS of Eq. (|A1[) . 

By analysing the minima of the density functional one can scan the free energy landscape and find whether there 
exist metastable states or not. Such a plan has been indeed followed using a simple approximation for J-"({p(x)}), 
called the Ramakrishnan-Youssouf density functional. It consists in only retaining, in the diagrammatic part of 
JF({/9(x)}), the linear and quadratic terms in p(x) — pa where po is the average density of the system. It reads: 



Fry({p{-x)}) = Hpo) + -s / dxp(x) log(p(x)/p ) - o / dxdy(>(x) - p )c(x - y)O(y) - p ) 
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where S(k) 



where c(x — y) is the directed correlation of the liquid, defined in Fourier space as c(k) 

is the structure function of the liquid, i.e. S(k) = (p(k)p(-k)}. 

The functional Try has been studi ed in great detail. It allows one to develop a quantitative theory of crystallisation 
by restricting to periodic solutions [207| | and, concerning liquids at low temperatures, to reveal the existence of 
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metastable amorphous solutions, i.e. glassy metastable states. This route was pioneered by Wolynes and Stoessel 
[72| . This first work has been followed up recently by many others [73j]. It is an approach that has the advantage 
of being very direct and concrete. Its drawbacks are that quantitative computations are difficult and computing 
properties like the configurational entropy is daunting, both numerically and analytically. 



2. Boundary Pinning Field and Replicas 



In the following we shall describe an approach based on the replica method which, although more abstract has the 
clear advantage that analytical computations for finite dimensional systems, even realistic ones such as hard spheres 
or binary Lennard- Jones mixtures, becomes feasible (although quite involved) [3] ■ 

To understand how replicas come about in a model without disorder, let assume that the system is in a regime of 
temperatures where there are indeed many very long-lived metastable states and that the Gibbs-Boltzmann measure 
is distributed over all of them, with the weight corresponding to the corresponding free energy, like in mean field 
models (see section Hi B II) . In order to study the statistical property of a typical metastable state, we focus on a very 
large cavity of radius R, carved in an otherwise infinite (or much larger) system. As explained in intuitive terms in 
section HlI B[ the basic idea is to apply a suitable boundary external field in an attempt to pin the system in one of the 
possible metastable states 236}. Contrary to simple cases, e.g. the ferromagnetic transition for which a positive or 
negative magnetic field selects states, the external "field" (or analogously the boundary condition) one has to impose 
to select a given amorphous state, is as unpredictable as the state one wants to select. To overcome this difficulty 
one can take an equilibrium configuration a, freeze the position of all particles outside the cavity, and use this as a 
boundary condition [99(. If the system is a thermodynamic glass characterised by many metastable states then this 
boundary condition should force the system inside the cavity to be into the same metastable state as the equilibrium 
configuration a. An illustration of this technique to the simple case of the Ising model helps to clarify and make 
concrete this discussion. We invite the reader to look at Appendix B for this purpose. 

Concretely, the procedure consists in computing the cavity partition function, for a fixed a: 



Z a (R) = J2^M-l3H(C))S(q out (C,C a ) - I) 



(A3) 



where q out (C,C a ) is a suitably defined overlap that measures the similarity between density configuration C and that 
of the a state in the space outside the cavity. When the overlap equals unity the two configurations are the same 
outside the cavity. In the large R limit, the intensive free energy of the metastable state, obtained by taking the 
logarithm of the partition function, is expected to be self-averaging and independent of C a . Physically, this means that 
the overwhelming majority of the metastable states sampled by the equilibrium Boltzmann measure are characterised 
by the same intensive free energy. 

Although we started from a problem without quenched disorder, we find that the analysis of the metastable states 
leads us to a problem where the configuration C a plays the role of a (self-induced) quenched disorder. In order to 
proceed further and compute the intensive free energy of a typical metastable states we have therefore to average over 
C a . As usual for quenched disorder problem, one can make use of the replica trick: 
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In order to compute the average (Z™ ) a one can introduce replicated configurations and write: 
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As usual with replicas one computes the above sum for integer and positive values of to — 1 and then makes an 
analytical continuation to make the m — * 1 limit. It is important to notice that the average in the numerator of the 
previous expression can be rewritten as the partition function of m replicas constrained to be identical outside the 
cavity but free to fluctuate inside since in the above expression C a is no longer different from the other replicas. 

Let us denote the logarithm of the partition function of the to constrained replicas as —(3F m . Once this quantity 
is known, one can compute the partition function of the large cavity as: 



J- = —T(\nZ a ) a = lim 



(F m — Fi] 
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dF m 
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This gives the free energy of one typical state inside the cavity. Since we are interested in thermodynamic quantities, 
henceforth we will consider R to be very large. If there are many states, i.e. an exponential number in the size of the 
system, then the free energy of the cavity T may be different from the unconstrained one F\ . This can be seen by 
rewriting the replicated partition function as a sum over all states with their Boltzmann weight. If the constraint is 



strong enough to force the replicas a — 1, . . . , m — 1 to fall into the same state as a itself, then: [237 



F m lnE„e-W° rajV In f df exp(N[-0fm + a(f,T)}) „ . x lk N 

— = -T N = -T — - — - — v 1 h J — u n/ = f m — TS(T), (A7) 

where /* is the free energy density that maximises the argument of the exponential, S(T) = a(f*,T) and N is the 
number of particles inside the cavity. Using Eq. (|A6[) . one immediately finds that /* = J-/N is the intensive free 
energy of a typical metastable state. But the free energy of the system without constraint is F\ = f* — T£(T) that 
contains the configurational entropy contribution. The replica method allows one to obtain both quantities, which 
then yields the configurational entropy: 



£(T) = -/3F I + ( 3^= 
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Hence, we have found that computing the statistical properties of metastable states reduces to the computation of 
the thermodynamics of m — > 1 replicas with the constraint that the overlap outside a spherical cavity of radius R is 
equal to one. In practice one has to do a computation for m replicas and take the space dependent overlap q a ,b{f) 
between them in the bulk of the cavity as an order parameter. Hence, one has to compute as accurately as possible 
the free energy as a function of the overlap q a ,b{r), and then find the stationary points. 

One always finds a trivial solution with uncoupled replicas, i.e. q a ,b(r) — for a ^ b. This is expected since if it 
was not for the boundary condition the replicas would be indeed completely uncoupled. If this is the only solution 
then F m = rnFi, and S = as it should be. One has therefore to inspect whether another solution exists. If it is the 
case then the constraint outside the cavity plays the role of a boundary condition for g a ,b(f) and selects this coupled 
replica solution. 

The approach outlined above is the starting point for several investigations of the glass transition that explicitly 
encapsulates the tenets of RFOT. It provides a microscopic basis to the arguments of section UlIBI on entropy driven 
cavity melting leading to the mosaic state. We have explained that one has to compute as accurately as possible the 
m-replica free energy as a function of the overlap q a j, and to look for a coupled replica solution. When the system is 
translation invariant one naturally looks for homogeneous solutions and this is indeed what has been done to compute 
quan titatively the configurational entropy and the properties of the metastable states and of the gla ss phase, see e.g. 
|208j | and [lj] for a review. But the study of inhomogeneous solutions is also very interesting. In |209l . |2 10j ] it has 
been shown that between Td and Tk (as computed by using the homogeneous solution), one finds an inhomogeneous 
solution with an overlap equal to one on the boundary that vanishes away from it on distances greater than a certain 
£*. This shows that the coupled replica solution actually is unstable and the metastable order does not propagate 
further than the length &* , which can be therefore identified with the mosaic length-scale discussed heuristically in 
section fill Bl This provides a microscopic way to compute the exponent 6. Up to now, all comput ations ha ve lead to 
the = d — 1, which is different from the value d/2 coming from KTW's wetting argument, see [I209j. |21(| for more 
details. 



3. Application: the "Weak Glass" expansion and the nature of the glass transition 

In order to apply the procedure described in the previous section to a real glass-former, one has to be able compute 
the thermodynamic properties of m replicas with the condition that their overlap is one outside a very large region 
S of linear size R. This leads to study a multi-component system characterised by an infinitely strong attraction 
between replicas outside S and no inter-replica interaction inside S. Same replica p articl es interact everywhere via 
the physical potential V. Using the diagrammatic formalism of liquid state theory [21 lj ] for this multi-component 
system one can write the Helmholtz free energy, which we will denote with a slight abuse of notations by F m , as a 
functional of the densities p a (x) and correlation functions C ab (x, x') = S(x — x") J2j <5( x ' ~~ x j))- As usual, it is 
more practical to use the following representation for the correlation function: 

C ab (x, x') =pl{\ + h a \x, x')) + PQ 5 ab 5(x - x') (A9) 

where po is the density of the liquid we are focusing on. The general expression of F m is quite complicated [2 111 . |212| | . 
In our case it can be simplified using explicitly that there is no interaction between replica inside S and assuming 
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translation invariance. This last assumption is certainly true far from the boundaries and simplifies the discussion 
without losing too much generality. 

As discussed above, the absence of metastable states translates into a solution where the replicas are completely 
uncoupled: = Po an d h ab — h8 a b- Instead, the coupled solution corresponds to = po but h ab = hi8 a b + h g (l — S a b) 
and allows one to unveil the existence of metastable states. In the m — * 1 limit, hi corresponds to the liquid correlation 
function and h g to the so-called called non-crgodicity factor q* : the plateau value of the density-density time correlation 
function, which measures the fraction of frozen-in density fluctuations. Note that contrarily to the case of spin glasses 
where replica symmetry breaking is spontaneous [2|], the boundary conditions break explicitly this symmetry in the 
present case. 

In the large R limit we find: 



where p is the chemical potential, the hat on h ab (k) denotes the Fourier transform, A the thermal Dc-Broglie wave- 
length, L(u) the function — log(l + u) + u — u 2 /2 and T the sum of all the more-than-doubly connected diagrams 
which are composed of po circles and h ab lines. 

It is clear that the replica uncoupled ansatz leads to a RHS which is not a function of m and, hence, to a zero 
configurational entropy via Eq. l|A8p . The replica coupled ansatz instead leads to a non zero configurational entropy. 
However, in order to obtain the correct result one has to minimise F m with respect to po and h ab , plug back the 
solution into the expression of F m and finally take the derivative with respect to m, which means performing at some 
previous stage the analytic continuation in m. Of course this is im pos sible to perform exactly and one must resort to 
approximations. This is what has been done in [20811212] ] . see also |l4j . In the following we shall present what we call 
a "weak glass" expansion with the purpose of showing that the glass transition is necessarily first order in h g despite 
the fact that it is second order thermodynamically. 

The weak glass expansion assumes that the boundary pinning field induces only a very small off-diagonal correlation, 
i.e. that the frozen-in density fluctuations h g are very small. We do not discuss the physical conditions that make 
this assumption self-consistent but surmise that it is possible to find a (presumably long-ranged) potential V such 
that this is the case. One can collect all terms which are of lowest order in h g keeping at the same time the hi 
contributions to all order. The zeroth order in h g is by construction independent of m. It coincides with the free 
energy for the one-component liquid we are focusing on and corresponds to the usual free energy of diagrammatic 
liquid theory as a function of hi. Analysing the expansion of T in h g , one finds that the first contribution containing 
h g is necessarily at least third order. Since all the other terms contain linear and quadratic terms in h g , we can drop 
the T contribution, which is sub-leading within the weak glass expansion. The other important observation, that 
can be checked straightforwardly, is that the terms containing h g are 0[m — 1) in the limit m — > 1. This has two 
important consequences. First, one recovers the usual liquid theory diagrammatics for m — ► 1, as one should. Second, 
it reveals that a non zero h g lead to a non trivial dependence on m of the RHS of (|A10[) and, hence, to a non zero 
configurational entropy via Eq. (|A8|) . 

As a summary, we find that the minimisation of F m in the limit m — > 1 leads to two independent variational 
problems. First, the usual one of liquid theory, from which one obtains the exact solution for the liquid h(x) and po. 
The second one gives access to the properties of metastable states in which the liquid can remain trapped for a long 
time. It corresponds to the minimisation of the m — 1 part of F rn . This is a functional of hi, h g ,pQ. One has therefore 
to plug the values of hi and po obtained from the solution of the first variational problem and minimise the result 
with respect to h g . The minimisation equation determining h g reads, to lowest order: 



where S(k) is the exact structure factor of the liquid. This shows that a continuous transition leading to a non zero 
h g is impossible except if the structure factor diverges for some vector (s) k, but this is not observed in experiments, 
and would correspond to a quite different physical situation. 

Therefore, one concludes that when a non zero h g (k) appears it necessarily does so in a discontinuous fashion. 
Thus, we find one of the very remarkable properties of glass-forming liquids which distinguish them from spin-glasses: 
when the dynamics starts to slow down above the glass transition, i.e. when metastable states emerge, the plateau 
in density-density correlation function appears discontinuously. This happens at a temperature which is identified 
with the dynamical temperature T c i as discussed in sect III B 31 . In order to understand quantitative the behaviour of 





46 



h g as a function of temperature one should go beyond the weak glass expansion and re-sum in some way an infinite 
subclass of diagrams. It would be interesting to pursue further this weak glass expansion. This approach would be 
the counterpart of density functional theory for crystallisation where one makes an expansion of the exact functional 
in powers of the p(x) — po- It may shed light on many important aspects of the dynamical transition that are very 
poorly known in finite dimensions. 

APPENDIX B: A TOY MODEL OF ENTROPY DRIVEN CAVITY MELTING 

In the following we shall focus on a toy model of entropy driven cavity melting. The aim of this section is to 
make more concrete several ideas explained previously and to shed some new light on procedures and intuitive 
ideas introduced in the study of the mosaic state. The model we focus on is the one dimensional Ising spin glass 
model with quenched random couplings J;, equal to +J and — J with probabilities 1/2. A gauge transformation 
Si — > TTj = o{Ji/ J)Si makes the model equivalent to the pure ferromagnetic Ising model. Thus, at zero temperature 
it is characterised by two degenerate ground states. The zero temperature order does not survive to any finite 
temperature. The mechanism that destroys the order provides a cartoon version of the entropy driven cavity melting 
discussed in section ITlI Bl In the following we will discuss this model quite in detail from this perspective. 

1. A recap on the low temperature physics of the one dimensional spin glass Ising model 

The Hamiltonian of the one dimensional Ising model reads: 

H = —^^JiSiSi+i Si = ±1 

i 

At T — and with open boundary conditions the model presents two degenerate ground states corresponding re- 
spectively to all spins up and all spins down after the gauge transformation. Reversing all the spins on an interval of 
size I costs an energy 4 J. This is independent of the length £ of the interval. In fact one can interpret this process 
as the creation of two defects. These are called kinks or domain wall and each one of them separates a region with 
all spins up from one with all spins down. One can considers the low temperature phase of the ID Ising model as a 
dilute gas formed by these excitations. Since each excitation costs an energy 2 J the density of kinks is p oc e~ 2J / T . 
This implies that a typical equilibrium configuration of the spin glass Ising model at low temperature consists - after 
gauge transformation - of up and down regions which alternate along the ID lattice on the scale Lk — p , see Fig 
1101 Note, however, that in the original gauge, the spin configurations look completely amorphous and without any 
sign of growing long range order. The spatial average of the connected correlation functions between two spins at 
distance x reads: 

i,j;j—i=% i,j;\i—j\=x fc=» 

In the large N limit it becomes zero except for x = 0, i.e. the configuration seems to lack any kind of long range 
order. Moreover, the Metropolis dynamics of this system become slow at zero temperature since it is governed by the 
diffusion of these rare defects. For instance, the time between two consecutive flips of a spin is of the order of the 
square of the distance between defects, i.e. T e 4J / T where r is the time on which each spin attempts to flip with the 
Monte Carlo dynamics. 

2. A mosaic with few tiles 

The low temperature configuration of the ID Ising model shown in Fig. fTU] can be interpreted as a mosaic state 
composed by just few tiles. There is indeed a strong similarity with the glassy mosaic. Finding back the physics 
described above from the mosaic point of view is an interesting exercise as we shall show in the following. 

Let us follow the procedure outlined for RFOT in section UlI Bl For simplicity we will work always directly on the 
pure Ising model, i.e. after having performed the gauge transformation. The mean field analysis of the problem is 
often related to the Curie- Weiss mean field equations: m = tanh 2/3 Jm, which would predict a phas e transition at 
temperature T c = 2 J. However, one can also write down the equivalent of the TAP free energy [213j . as one would 
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FIG. 10: Top: typical low temperature configuration for the Ising spin glass at low temperature. Bottom: corresponding 
configuration after gauge transformation. Lk is the typical distance between kinks. 

naturally do in a RFO T analysis. In this ID case the free energy as a function of the local magnetization can be 
computed exactly [2131 ] - For simplicity we just reproduce its continuum approximation which reads: 

F[{m(x)}] = I dx [D{X7m{x)) 2 + g(m 2 {x) - l) 2 ] 



where D and g are two microscopic constants. The analysis of the low free energy minima reveals that there are two 
absolute minima corresponding to uniform plus or minus magnetisations. Moreover, there are many other minima 
that corresponds to have L domain walls inside the system. Their free energy (or energy since they coincide in the 
low temperature limit we are focusing on) is higher than the ground state one of an amount AF = 2JL. One can 
easily compute the number of such states: it corresponds to the different ways to place the kinks on the lattice and 
reads n^jgwj • This l ea ds to an intensive 'configurational entropy' equal to S c — — plogp — (1 — p) log(l — p), where 
p = L/N is the density of the kinks. Note that for the Ising case, the minima are not separated by barriers contrary 
to what one expects for the glassy minima of super-cooled liquids. 

Following the same procedure described in the previous section for glasses, we compute the partition function of the 
system by summing over all states and weighting each term by its Boltzmann weight: 

a 

where Fqs is the free energy for the state without kinks. The integral can be performed by the saddle point method. 
The value of p that dominates the integral at low temperature is given by p* ~ exp (— %r). So we find that the 
'configurational entropy', S c (p*), is small and goes to zero at T = 0. Note, however, that real space configurations 
still looks completely disordered. Only applying the gauge transformation one can discover a growing order. Thus, 
at low temperature we are in a situation very similar to glasses: the spatial average of the two point correlation 
function is zero except at zero distance, indicating apparently no growing order, and at the same time the dynamics 
slows down very fast. Using the argument developed in section fill Bl we would expect that the competition between 



configurational entropy and surface tension leads to a mosaic state with a length (j^rj • I 11 this simple ID 

case d = 1,6 = d — 1 = and it is natural to take T = 2 J (or more generally a constant times J). This leads to an 
estimate of the mosaic length e^F . 

This result can substantiated by computing explicitly the average overlap inside a frozen cavity, as described in 
section [Til Dl and Appendix A. Let's consider an equilibrated configuration {S^ 9 } and freeze all the spins at a distance 
larger than R from a given spin Si- We analyse the thermodynamics of the cavity with this boundary condition, in 
particular we compute the average value of Si. It is easy to check that this is given by: 

(s) = gVj nUw tanh/3Jfc + gkHS^ tanh/3Jfc 
{sr} ~ i+^ +4 sr%nt + r^tanh/3j fe 



l/(d-9) 



As expected, this average value depends on the equilibrated configuration only through its boundary spins. As in 
Appendix A, we shall now compute the average overlap between (5 l i){s E 9} and the value of the spin S^ q in the original 
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equilibrated configuration: 

1 {S° q } 

Note that we also average over the site i. The results will be therefore equivalent to averaging over the quenched 
disorder. The computation of q(R) is rather straightforward and will not be detailed here. The final result is: 

2(tanh/3J) 2 « 
I[ ' l + (tanh/3J) 2 « 

The function q(R) decreases monotonically on a length-scale £ = 2 intanh ) 3j - ^ n tne ^ ow temperature limit we are 
interested £ ~ e 2,///T /4. This coincides with Lk up to a constant factor. Thus we find the same results than in the 
previous heuristic argument and we conclude that the system is indeed in a mosaic state of length Lk- Furthermore, 
we find that the configurational entropy of the system is equal to S c (p*). Since the intrastate entropy is zero this is 
also the entropy of the system and indeed it coincides with the correct result for the ID Ising spin glass. 



3. Further analogies and differences 

We have found that the procedure outlined in Appendix A to study the mosaic state successfully applies to the toy 
model considered here. Of course there are important differences between the ID Ising spin glass and the behaviour 
conjectured by RFOT for a glass-forming liquid. The most important ones are that (i) metastable states are not 
separated by barriers that grow with size, actually they are connected by zero modes; (ii) the configurational entropy 
on the mosaic length scale £ grows only logarithmically with £, contrary to three dimensional glasses where it is 
expected to scale as £ e with 3/2 < 9 < 2. This means that the mosaic of the Ising spin glass has much less tiles than 
the one conjectured for 3D supercooled liquids, (iii) Although linear responses and two point correlation functions do 
not show any long range order, three point responses and the square of the connected correlation function do so, as 
in spin glasses. In glasses this is expected to take place at finite frequency only. 

Despite these differences the analysis of the toy model is interesting because it allows one to make concrete many 
concepts and procedures. It could also be used to test new ideas, as we do in the following. 

Clearly, RFOT is rooted in the analysis of mean field glassy systems. However, as we have seen, its finite dimensional 
extension departs in very important ways. In order to develop a complete and self-consistent theory, it seems necessary 
to develop an approach which does not make too much reference to mean field systems. In many instances, there are 
observables and concepts that are well defined within the mean-field approximation only. Their extension to finite 
dimensional is somehow fuzzy and make the whole RFOT construction shaky. A very good candidate would be a 
renormalisation group description of the ideal glass transition. Let us use the toy model to have a first hint of what the 
outcome would be. In the ID spin glass Ising model the real space renormalisation group can be performed exactly. 
Using the Migdal-Kadanoff prescription, one integrates out the spins on the, say, odd sites to get a new effective 
coupling between even spins. The exact RG relation reads: 

eff 1 1 + tanh/jjj tanh/3J m 
P 4 ,*+2 - 2 n 1 -tanh/JJitanh/SJi+i 

By iterating the RG transformation one finds that 0J e ff becomes of the order of one on scales of the order Lk- This 
means that on this scale (and above) the system is effectively at high temperature: it has almost no stiffness and the 
entropy is of the order of one. This suggests an RG way to characterise the mosaic state in liquids. Integrating out 
degrees of freedom up to scale £, one should find that the entropy increases, the surface tension between metastable 
states decreases, and the relaxation timescale increases. When £ equals the point to set length one should find a 
liquid characterised by a vanishing surface tension between states. Thus, metastable states should become unstable 
and the related configurational entropy disappears. Clearly, the above scenario is highly speculative but we do hope 
that future works will be performed along these lines. 
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53 



corresponds in that case to the shear modulus for frequencies much larger than r" 1 , see section [ill E 
[222] The whole cavity scenario has been tested in details in the context of "usual" nucleation theory, see 
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and not between two typical TAP states, is a little different from ours above. 
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conclusions. 

[230] Here there is a subtlety related to the existence of conserved quantities in the dynamics. If this is the case, as it for 
real liquids, then the upper critical dimension is d c — 8. Instead, without any conserved quantities, as it is the case for 
the p-spin models with Langevin dynamics for example, the upper critical dimension is 6. The validity of the previous 
argument, which assume scaling, could therefore be questionable. Nonetheless a more careful treatment 96] shows that 
the final result is corre ct. 

[231] We refer the reader to [1291 . Il82l . Il83l . Il84l . Il85| for other interesting ideas about the glass transition. 

[232] It would in fact be interesting to measure Up in the same conditions. A breakdown of the relation between this observable 

and the relaxation time as R is varied may indicate that Eq. (|72[) is less fundamental than anticipated. 
[233] Note that for v — 2/3, this predicts a simple linear behaviour of S xa (T) below T* . 

[234] Alternative words for this are: "jammed", "immobile" or "inactive", but all these ideas reflect the fact that the relevant 

configurations are at least locally stable. 
[235] Google scholar counts 1882 citations of the Adam-Gibbs paper on November 5th, 2009, since the beginning of the nineties. 
[236] Our presentation is different from that of the original papers [541 . [78| but based on the same ideas 

[237] The alert reader should be very suspicious about this assumption, and rightly so. This will be critically reconsidered later 
in this Appendix. 




